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ABSTRACT 

We find a new mechanism of neutron star radiation wherein radiation is produced by the stellar 
interior. The main finding is that neutron star interior is transparent for coUisionless electron 
sound, the same way as it is transparent for neutrinos. In the presence of magnetic field the 
electron sound is coupled with electromagnetic radiation, such collective excitation is known as a 
fast magnetosonic wave. At high densities the wave reduces to the zero sound in electron liquid, 
while near the stellar surface it is similar to electromagnetic wave in a medium. We find that 
zero sound is generated by superfluid vortices in the stellar core. Thermally excited helical vortex 
waves produce fast magnetosonic waves in the stellar crust which propagate toward the surface 
and transform into outgoing electromagnetic radiation. The magnetosonic waves are partially 
absorbed in a thin layer below the surface. The absorption is highly anisotropic, it is smaller 
for waves propagating closer to the magnetic field direction. As a result, the vortex radiation is 
pulsed with the period of star rotation. The vortex radiation has the spectral index a k, —0.45 
and can explain nonthermal radiation of middle-aged pulsars observed in the infrared, optical and 
hard X-ray bands. The radiation is produced in the star interior, rather then in magnetosphere, 
which allows direct determination of the core temperature. Comparing the theory with available 
spectra observations we find that the core temperature of the Vela pulsar is T w 8 x lO^K, 
while the core temperature of PSR B0656-I-14 and Geminga exceeds 2 x lO^K. This is the first 
measurement of the temperature of a neutron star core. The temperature estimate rules out 
equation of states incorporating Bose condensations of pions or kaons and quark matter in these 
objects. The estimate also allows us to determine the critical temperature of triplet neutron 
superfiuidity in the Vela core Tc = (7.5 ± 1.5) x lO^K which agrees well with recent data on 
behavior of nucleon interactions at high energies. We also find that in the middle aged neutron 
stars the vortex radiation, rather then thermal conductivity, is the main mechanism of heat 
transfer from the stellar core to the surface. The core radiation opens a possibility to study 
composition of neutron star crust by detection absorption lines corresponding to the low energy 
excitations of crust nuclei. Bottom layers of the crust may contain exotic nuclei with the mass 
number up to 600 and the core radiation creates a perspective to study their properties. In 
principle, zero sound can also be emitted by other mechanisms, rather than vortices. In this case 
the spectrum of stellar radiation would contain features corresponding to such processes. As a 
result, zero sound opens a perspective of direct spectroscopic study of superdense matter in the 
neutron star interiors. 



1. Introduction 



Properties of matter at densities much larger than nuclear are poorly known and constitute a challenging 
problem of modern science. They have broad implications of great importance for cosmology, the early 
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universe, its evolution, for compact stars and for laboratory physics of high-energy nuclear collisions. Present 
searches of matter properties at high density take place in several arenas. One of them is investigation of 
neutron stars. 

According to theory of neutron stars (NSs) (Yakovlev et al. 1999a) a NS can be subdivided into the 
atmosphere and four internal regions: the outer crust, the inner crust, the outer core and the inner core as 
shown in Fig. 1. The atmosphere is a thin plasma layer, where thermal electromagnetic radiation is formed. 
The depth of the atmosphere varies from some ten centimeters in a hot NS down to some millimeters in a cold 
one. The outer crust extends from the bottom of the atmosphere to the layer of density ~ 4.3 x lO^^g/cm^ 
and has a depth of several hundred meters. Ions and electrons compose its matter. The depth of the inner 
crust may be as large as several kilometers. The inner crust consists of electrons, free neutrons and neutron- 
rich atomic nuclei. Neutrons in the inner crust are in a superfluid state. The outer core occupies the density 
range 0.86 — 2^0^ where tpo = 2.8 x lO^^g/cm^ is the saturation nuclear matter density, and can be several 
kilometers in depth. It consists of neutrons n with some (several per cent by particle number) admixture 
of protons p and electrons e. Almost all theories of NSs predict the appearance of neutron superfluidity 
and proton superconductivity in the outer NS core. In a low-mass NS, the outer core extends to the stellar 
center. 

More dense NSs also possess an inner core (see Fig. 1). Its radius may reach several kilometers and 
central density may be as high as 10 — 15ipo- The composition and equation of state of the inner core are 

poorly known and it constitutes the main NS "mystery" . Several hypotheses are discussed in the literature. 
One of them is appearance of S- and A-hyperons in the inner core. The second hypothesis assumes the 
appearance of pion condensation. The third hypothesis predicts a phase transition to strange quark matter 
composed of almost free u, d and s quarks with a small admixture of electrons. Several authors also considered 
the hypothesis of kaon condensation in the dense matter. The radius of the NS usually decreases as the 
stellar mass increases (NS becomes more compact). The mass reaches its maximum at some density, which 
corresponds to the most compact stable stellar configuration. Typically the stellar mass ranges between 
0.5 — 3Mq for different equations of state {Mq is the solar mass). 

Despite of 30 years enormous theoretical and observational efforts, the properties of matter in neutron 
star cores remain poorly known. The main reason is that detected thermal radiation from the star surface 
carries practically no information about the stellar interior. Some indirect methods, such as estimates 
of neutron star masses, radii and cooling rates for different equation of states, have been applied to find 
constrains on the matter behavior. However, non of them is reliable enough because the answer strongly 
depends on the model parameters which are also poorly known. 

In this paper we discuss a new direct way to determine properties of the superdense matter. We found 
that neutron star interior is transparent for electron zero sound, the same way as it is transparent for 
neutrinos. Zero sound is a collective density (or spin-density) wave in Fermi liquids that was predicted by 
Landau in 1957 and later observed in ^He and electron liquid in metals. The wave frequency w is greater 
then the frequency of interparticle collisions which makes the wave different from usual hydrodynamic sound. 
In neutron stars the zero sound can propagate in electron fluid which flUs the stellar crust and core, the 
collisionless (zero sound) regime is reached at a;/27r > lO^^Hz. The large density of electrons in NS interior 
results in substantial reduction of the interaction energy between fermions as compared to their kinetic 
(Fermi) energy, the electron gas becomes more ideal. In such weakly interacting Fermi system the attenuation 
of zero sound is exponentially small (see Eq. (A9) in Appendix A). As a result, zero sound propagates across 
the star in a ballistic regime, which is similar to neutrino. However, this is not the case for single-particle 
excitations in electron liquid, their mean free path usually does not exceed 0.01 — 1cm. The single-particle 
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excitations are in a local thermodynamic equilibrium with the stellar matter, while collective zero sound 
modes do not. At the same time, negligible attenuation of the zero sound creates an obstacle to excite such 
waves inside the stellar medium. Some special processes are required for their generation. We found that 
zero sound is generated at the core-crust interface by helical motion of superfluid vortices. Vortex modes 
themselves arc excited thermally via single-particle channel, they are in thermodynamic equilibrium with 
the surrounding medium. 

In the presence of magnetic field the electron sound is coupled with electromagnetic radiation, such 
excitation is known as a fast magnetosonic wave. At large densities the effect of magnetic field is negligible 
and the wave behaves as a sound wave, however near the stellar surface (small matter density) it reduces to 
electromagnetic wave which leave the star and propagates into surrounding space. Such property results in 
a new mechanism of star radiation produced at the stellar interior. Contrary to neutrino, electromagnetic 
radiation is easy to detect which opens a perspective of direct spectroscopic study of superdense matter. 
Thermally excited helical vortex waves generate fast magnetosonic waves in the stellar crust which propagate 
toward the surface and transform into outgoing electromagnetic radiation. We found that such mechanism 
results in radiation with non Planck's spectrum (Eqs. (23), (42)) and theoretical predictions agree well 
with nonthermal radiation of some middle-aged pulsars observed in the infrared, optical and hard X-ray 
bands (Figs. 5, 6). The important result of our theory is that detection of vortex radiation in the X— ray 
band allows direct determination of the core temperature. From available spectral data we found that the 
core temperature of the Vela pulsar is T « 8 x lO^K, while the core temperature of PSR B0656+14 and 
Geminga exceeds 2 x lO^K. This is the first measurement of the temperature of a NS core. The cooling 
rate and temperature of the central part of a NS substantially depends on the properties of dense matter. 
Our temperature estimate excludes exotic equation of states incorporating Bose condensations of pions or 
kaons and quark matter in these objects. Such conclusion agrees with recent measurement of gravitational 
redshift of a neutron star (Cottam et al. 2002). Moreover, the temperature estimate allows us to determine 
the temperature of superfluid transition of neutrons in the Vela core Tc = (7.5 ± 1.5) x lO^K. This value 
agrees well with theoretical predictions that are based on new nucleon-nucleon interaction potentials which 
fit the recent world data for pp and np scattering up to 500Mev (Baldo et al. 1998; Arndt et al. 1997). 
Examples of such potentials are Nijmegen I and II. For the matter density 2'tpo (which is a central region 
density of a canonical NS without exotic inner core) these nuclcon potentials predict Tc ~ 7.5 x lO^K for 
neutron superfluidity with (triplet) pairing (Baldo et al. 1998). One should mention that in order to 
find the transition temperature at densities up to lOt/'o more accurate potentials which fit the measured 
nucleon-nucleon phase shifts up to about IGev are required. A systematic study of pulsar radiation in the 
V— ray band can locate objects with fast cooling core which could be candidates for NSs with exotic states 
of matter. 

The plan of the paper is the following. In section 2 we discuss neutron vortices in the superfluid core 
and generation of magnetosonic waves by vortices at the crust-core interface. In sections 3, 4 we investigate 
propagation of magnetosonic waves across the stellar crust and their transformation into electromagnetic 
radiation at the surface. In section 5 we consider spectrum of electromagnetic radiation produced by vortices 
and analyze the applicability of our results. In section 6 we compare the theory with radiation spectrum 
of some middle-aged pulsars observed in infrared— X-ray bands and discuss properties of superdense matter 
based on the observed spectra. Section 7 is devoted to a problem of heat transport across the stellar crust. 
In section 8 we draw our conclusions and discuss perspectives in the field. In Appendices A, B, C we consider 
zero sound in an electron liquid, derive equations of coUisionless magnetic hydrodynamics for electron motion 
in the stellar crust and investigate propagation of electromagnetic waves across the atmosphere. 
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2. Vortex radiation in rotating neutron stars 

If a NS rotates with an angular velocity CI =ftz a periodic vortex lattice forms in a superfluid neu- 
tron phase of the stellar core. Typical value of the lower critical angular velocity of the vortex formation is 
Oci/27r ^ 10^^"' .s^^ (Sedrakyan & Shakhabasyan 1991). The vortices are parallel to the axis of rotation apart 
from a tiny region (< 10~'^cm) near the condensate surface, where they are curved and cross the boundary 
in the perpendicular direction. The distance between vortices in a triangular lattice isb= {Airh/VSMfiy/'^, 
where M is the mass of the condensate particles. E.g., for a neutron superfluid (M — 2m„, where m„ 
is the neutron mass) and Q/2n = 10s~^ the lattice period 6 = 6 x 10~^cm. For NSs the characteristic 
size of the vortex core is ^ = 10~^^cm, so that b ~ 10^°^ and the total number of vortices in the super- 
fluid A'' = MR'^Cl/U Ri 10^'', where R ~ 10km is the radius of the superfluid condensate (Sedrakyan & 
Shakhabasyan 1991). 

First, let us consider a single straight vortex line with length L located in a spherical one component 
superfluid condensate. We assume the condensate density is uniform inside the sphere. One should note 
that this situation differs from real NS core which consists of a mixture of neutrons, protons and electrons. 
However, main properties of vortices can be understood within the simplified model. Normal modes of the 
vortex are helical waves (Lifshitz & Pitaevskii 1980) 

r = Aexp{iu;t + ik\iz), a; = ^ A;| ln(l/A;||^), (1) 

where A is the amplitude and fc|| is the wave number. Eqs. (1) are valid for <C 1, or w <C Wc = Ti/2M^^ w 
lO^^Hz and describe rotation of a helix with the angular frequency uj. 

To find the amplitude of the helical waves A one can use the expression for the energy of the system in 
terms of normal vortex modes (Svidzinsky & Fetter 2000a). If a superfluid contains a vortex line of length 
L it changes the kinetic energy of the fluid by 

E, = JMl^d^r = nLn^Hb/a (2) 

where n is the superfluid particle density. Change in the kinetic energy (2) gives the main contribution 
to the change in the fluid energy. One should mention that the superfluid density is suppressed near the 
vortex core. However, the normal component fills the vortex core so that the total fiuid density remains 
approximately constant (Nygaard et al. 2002). This occurs because any variation in the total density in 
a Fermi liquid (under NS conditions) produces substantial change of the Fermi energy and increases the 
energy of the system. The same situation takes place during the vortex motion and, as a result, the helical 
vortex waves are not accompanied by oscillation in the total density. In the wave the superfluid and normal 
components oscillate in the opposite phase so that the total fluid density remains approximately constant. 

Taking into account that excitation of the helical wave (1) changes the vortex length by the value 
iA^fc|/2, one can obtain the energy of the system in terms of normal mode amplitudes (Svidzinsky & 
Fetter 2000a) 

E = Eo + TmLY^ Al^huJrn- (3) 

m 

In the NS core the helical vortex waves are excited thermally by particle scattering on the vortex cores. 
Vortex excitations are in a local thermal equilibrium with the stellar matter. Then at temperature T the 
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energy of the wave with the frequency w is equal to hui times the mean number l/[eyi\i{1iLo /UbT) — 1] of 
excitations in that state. Therefore, the vortex wave amphtude is given by (cf. (Barenghi 1996)) 



TrnL[exp{hu} / kBT) — 1] ' 

where fcs is Boltzmann's constant. For typical parameters of NSs A <^ ^. Indeed, if we take n = lO'^^cm"'^, 
L = 10km, a;/27r = lO^^Hz, T = lO^K and ^ = 10~^^cm, we obtain for the amplitude of the thermally 
excited helical wave A w 10~^^. 

Apart from helical waves there are excitations of the vortex lattice in a plane perpendicular to the 

vortex lines which correspond to relative displacements of vortices with respect to each other (the so called 
Tkachenko's modes). These excitations are analogous to a vortex precession in trapped Bose-Einstein con- 
densates where a vortex line moves along a trajectory with constant trapping potential (Svidzinsky & Fet- 
ter 2000b). However, frequencies of the vortex oscillation due to excitation of such modes are small. Indeed, 
the frequency of Tkachenko's modes is given by lot = \/fiQ/4:Mk±, where k± is the wavenumber in the 
plane perpendicular to vortices (Tkachenko 1965). For maximum frequency k± ~ 1/6, where b is the lattice 
spacing, and we obtain 

max ^5 which is much less than frequencies of radiation we are interested in. 



2.1. Energy transport by helical vortex waves in superfluids 

Let us consider a superfluid with a vortex line located along the z axis and introduce perturbation in the 
particle density n' and the velocity potential x'- i^ = iT'O + n' , x = Xo + x' t where no = no{r) is the superfluid 
density for a straight vortex line and xo = fi^/M. At distances r ^ (we use cylindrical coordinates r, (p, z 
measured from the center of the vortex line) the helical vortex motion (1) can be described in terms of n' 
and x' as 

n' = — ^exp [i{u;t + k\\z — (p)] (5) 

X' = iA^exp[i{u;t + kiiz-<f>)] . (6) 

In the bulk the function x' (as well as n') should satisiy a linearized equation of superfluid hydrodynamics, 
which is, in fact, the wave equation for sound propagation 

^ - clV\' = 0, (7) 

where Cg is the speed of sound in superfluid. In the limit u> <^ Uc one can omit the flrst term in Eq. (7). 
Then solution of this equation which satisfies the boundary condition (6) has the form 

X' = iA^K^ik^r) exp [i{ijt + k^^z - </))] , (8) 

where Ki is the modified Bessel function. The solution exponentially decreases at r » Ki{k\ir) ~ 

(7r/2fc||r)^/^ exp(— fc||r). We obtain that helical vortex waves produce perturbation in the superfluid velocity 
V = Vx which is localized near the vortex core. Such perturbation results in the energy flow along the 
vortex line so that the energy flux is localized near the vortex. 



The energy flux can be estimated as 



Q = llpV"^.dS, (9) 
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where p = Mn is the superfluid density, V' is the perturbation in the supcrfluid velocity due to hehcal vortex 
motion and doj/d'k =?ifc|| ln(l/A:||^)i/M is the group velocity of the helical wave. In Eq. (9) the integral is 
taken in a plane perpendicular to the vortex line. Then, we note that L J ^pV'^dS is the total energy of the 
helical wave, which at temperature T is given by 



pV'^dS = ^ (10) 



Taking Eq. (10) into account, we find the following expression for the energy flux along the vortex: 

LVM [exp{nu/kBT) - 1] ■ 



The number of vortex modes within the interval dk\\ is Ldk\\/2'K = L\fM duj / 2tt \/ Tiuj ln(a;c/a») and, hence, 
the spectral power density of energy flux along a single vortex line is 

^^'^^ = 27r[exp(;ia;/fcBT)-l]- ^^^^ 

This expression is a universal function in the sense that it depends only on the temperature. At Ttw < ksT 
the energy spectrum is independent of frequency and P = kBT/2-jT. 

One should note that the energy flow produced by helical waves substantially differs from energy radi- 
ation in the case of a moving straight vortex line. A moving straight vortex in a compressible fluid emits 
sound waves in the direction perpendicular to the vortex line. A homogeneous two-dimensional (2D) super- 
fluid is equivalent to (2-|-l)-dimensional electrodynamics, with vortices playing the role of charges and sound 
corresponding to electromagnetic radiation (Ambegaokar et al. 1980; Arovas & Freire 1997; Ovchinnikov & 
Sigal 1998; Fischer 1999; Lundh & Ao 2000). Thus, radiation of soimd waves from a straight vortex moving 
on a circular trajectory is concentrated in a plane perpendicular to the vortex line and can be estimated 
using formulas for the cyclotron radiation of an electrical charge moving along a circular orbit in 2D space. 

To understand the difference between energy radiation by a straight vortex and a helical wave we can 
solve the problem in general case. Indeed, general solution of Eq. (7) that satisfies the asymptotic behavior 
(6) is , 

x' = — ^j^k±Yi{k±r) exp [i{ujt + k\\z - </>)] , (13) 



where Yi is the Bessel function and k± = ^w^/c^ — At large r the solution (13) has the asymptotic 

, i^/TThA^/k± /, 37r\ r., , ,^ 
X = ^My/¥ ^'"^ V ~^)^^^ + <P)\ , (14) 

, n dy' ^/nhLunA r — / 37r\ r., , 

" = = - v2Mc^v^ ^'^" r^" - tJ [^(^^ + ^11" - "^^J • ('^^ 

For LO > Csk\i {k± is a real number) the solution describes radiation of sound waves with a conical wavefront 
(the cone angle is arctan(fc_L/A;||)). The particular case = corresponds to a circular motion of a straight 
vortex and sound waves are radiated in the plane perpendicular to the vortex line. However, in the case of 
helical waves the wavelength of sound that would be emitted is much larger than the wavelength of helical 
wave and the condition of sound radiation w > Cgfen (which for a helix is equivalent to a; > 4a;c/ln(a;c/w)) 
fails. As a result, the solution (13) reduces to (8) which describes energy flow along the vortex. 
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2.2. Magnetic field associated with neutron vortices 

Now let us take into account a superconducting proton component in NS core. The bulk of supercon- 
ducting protons do not rotate by forming vortices, which differs from neutron superfluid. Superconducting 
protons corotate with the crust and electron fluid by generating a surface current which produces a small 
uniform magnetic field in the interior H = —2mpcd/e, here mp, e is the proton mass and charge (Alpar 
et al. 1984). The vector potential associated with this field A = — mpc(0 x r)/e yields the proton velocity 
which corresponds to rigid body rotation 

while the phase of the proton order parameter tp is constant inside the superconductor. However, due 
to interactions between protons and neutrons, neutron superfluid velocity generates a superfluid current of 
protons in the vicinity of neutron vortices (drag effect) and an associated magnetic field around each neutron 
vortex line. The magnetic field of a neutron vortex is confined within a radius of the order of the effective 
London length = mpm*c^/47re^/9p, where pp is the mass density of superfluid protons and m* is the 
proton effective mass (Alpar et al. 1984). For the interior density p„ = 10^"'g/cm'^ and proton concentration 
Pp/ Pn = 0-05 the screening length is A^ 10~^^cm. At r > ^ the magnetic field is given by 

where = $o<5to*/to„, $o = i^fi-c/e is the flux quantum, (5m* = m* — nip is the contribution to the proton 
effective mass due to the interaction with the neutron medium and Kq is the modified Bessel function. 
Sjoberg (1976) has shown that for typical barion densities of a NS core 1 — 6 x lO-'^^g/cm^ the ratio m*/mp 
changes approximately from 0.5 to 0.25. This result is important for our theory since it indicates that 
neutron vortices are strongly coupled with protons that carry electric charge. As a consequence, motion of 
neutron vortices is accompanied by motion of protons and, therefore, electrons. This creates a possibility of 
excitation collective modes in the electron liquid by oscillating neutron vortices. In superconducting (npe) 
phase the magnetic field lines are parallel to neutron vortices which is similar to magnetic field of a long 
solenoid with the radius A^ . In the stellar crust there is no drag effect and magnetic lines diverge the same 
way as magnetic lines near the edge of the solenoid. 

One should note that fossil (not associated with neutron vortices) magnetic field can exist in the stellar 
core. However, the process of the magnetic fields generation in NSs and, hence, their structure and localiza- 
tion is still under discussion. There are arguments that NS magnetic fields may be confined to the crustal 
layer only, where it has been generated by thermomagnetic instability during the first years of NS evolution 
(Blandford et al. 1983; Urpin et al. 1986; Wiebicke & Geppert 1996). However, there is a possibility that 
magnetic field permeates both the stellar crust and the core. Such magnetic field (if it exists) forms before 
the transition into superconductivity by convection which arises in first 10-20 seconds after the NS is born 
and most probably, has a very complicated twisted structure (Thompson & Duncan 1993). As a rule the 
magnetic flux tubes inside a neutron star do not parallel to the axis of rotation (Ruderman et al. 1998). 
One should mention, that typical value of the critical flcld for vortex formation in superconducting core, 
lO^^Gs, is much greater than actual magnetic field in NSs lO^^Gs. Therefore, the core state with penetrating 
fossil magnetic field is metastable. Proton vortices are expelled from superconductor by different dissipation 
mechanisms which determine the life time of the metastable configuration. One of the mec;hanisms is the 
following. Since the neutron vortex and a fiux tube strongly interact as they pass through each other, the 
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moving neutron vortices will push on the proton's flux-tube array, forcing it to move outside in a spinning- 
down NS (Sauls 1989). As neutron vortex motion moves an entrained flux tube, that tube is ultimately 
pushed into the crust-core interface for almost any initial flux-tube conflguration. E.g., during the spinning- 
down stage when the initial stellar rotation period of 1ms (at the star birth) increases to 100ms (middle aged 
pulsars) the spacing between neutron vortices increases by a factor of 10. As a result all magnetic tubes 
should be pushed out into the crust unless originally some of them pass very close (at a distance less then 
0.1 stellar radii) to the axis of rotation. 

In this paper we will not consider possible presence of fossil magnetic tubes in the stellar core. Such tubes 
would not modify radiation produced by neutron vortices, however, they may give an additional contribution 
into stellar radiation. The main reason why wc do not consider the possible presence of fossil magnetic tubes 
is that the observed spectra of middle aged pulsars can be explained well only by radiation from neutron 
vortices. The other argument is that the core temperature we determined from the observed spectra is very 
high. The big difference between the core and surface temperatures suggests strong heat isolation of the 
stellar core. In Sec. 7 we show that geometry in which magnetic lines at the crust bottom are parallel to the 
core surface (or magnetic fleld has a closed configuration within the star) provides sufficient heat isolation 
and can explain the observations. Such geometry suggests that fossil magnetic field, if it has been generated 
in the core, already expelled from the superconducting core at the stage of middle aged pulsars. Based on the 
results of this paper we can provide another argument against presence of fossil magnetic field in the stellar 
core. The point is that our findings show that superfiuidity of neutrons and superconductivity of protons 
in the NS core are very strong. Such conclusion is also supported by the theory which is based on new 
realistic interaction potentials between nucleons (Baldo ct al. 1998; Elgaroy et al. 1996). The temperature 
of superconducting transition for protons near the stellar center is expected to be no less than 2 x lO^^K 
(Elgaroy et al. 1996). The modified Urea process cools the stellar core to such temperature during the time 
r Ri lyr/Tg = 0.5s (Pethick 1992). This means that stellar interior practically instantly after the star 
formation becomes superconducting and the convection mechanism cannot generate magnetic field at the 
stellar core simply because there is no time for that. 



2.3. Generation of magnetosonic waves at the boundary of superconducting phaise 

At the interface between the stellar crust and the outer core (at the density 2.4 x 10^*g/cm^) the 
matter undergoes a first order transition from a nucleus-electron-neutron {Aen) phase to a uniform liquid of 
neutrons, protons and electrons {npe phase). Nuclei in the Aen phase can be spherical and form a Coulomb 
crystal or have a deformed shape which results in an exotic distributions of nuclear matter. Transition 
from the Aen phase to the uniform liquid is accompanied by a small (about 1.4%) density jump (Douchin 
& Haensel 2000). In the npe phase the proton liquid is superconducting (Baym et al. 1969), while in the 
Aen phase protons constitute a part of nuclei. As a result, the Aen phase is not superconducting, however, 
neutron component is still superfluid. 

As we have shown before, helical vortex motion does not generate sound waves in the stellar core because 
the condition of sound radiation is not satisfled. However, at the interface between the superconducting (npe) 

and Aen phases such vortex motion produces magnetosonic waves which propagate across the stellar crust 
toward the surface. One should note that frequencies of radiation we are interested in arc larger than effective 
frequencies of electron-phonon and electron-electron collisions fe.ph, t'e.e- According to Eqs. (30), (45) at 
T = lO^K, p = 10"g/cm^ ne = 10^^-^ for ^^Fe lattice we obtain i^o.ph = 8 x 10^^s-\ ly^.c = lO^^s''^. So, 
at u/2Tr > lO^^Hz the electrons in the stellar crust behave as an independent system of Fermi particles and 
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move in collisionless regime. In such regime a zero sound is a possible collective excitation of Fermi liquid. 
In details the zero sound in electron liquid is discussed in Appendix A. In the presence of magnetic field the 
electron sound is coupled with electromagnetic radiation. In Appendix B we derive equations of electron 
motion in collisionless regime and show that they reduce to usual equations of magnetic hydrodynamics in 
which speed of zero sound enters the equations instead of speed of usual sound and the Alfvcn velocity is 
estimated in terms of the electron density pe = meUe'- ua = H/^JAirpe- One should mention that near the 
crust bottom the Alfven velocity is small because of large density. \{ H = lO^^Gs and pe = lO^^g/cm^ we 
obtain ua = 9 x lO^cm/ s which is much less than the speed of zero sound Us ~ c, where c is the speed of light. 
In magnetic hydrodynamics there arc three possible types of waves: Alfven, fast and slow magnetosonic. 
Here we consider radiation of fast magnetosonic waves which are in the limit Ug ^ ua reduce to zero sound. 
Vortex motion also generates Alfven and slow magnetosonic waves in the stellar crust. However, the intensity 
of radiation of those waves is small (in parameter u\/u^g) as compared to the intensity of fast magnetosonic 
wave. Also Alfven and slow magnetosonic waves are highly damped and can not reach the stellar surface. 
Due to these reasons we will not discuss radiation of Alfven and slow magnetosonic waves. 

Let us now consider excitation of electron zero sound by vortices in detail (see Fig. 2). Helical waves 
produce oscillation of the neutron vortex and, therefore, oscillation of protons in the npe phase. In Aen 
phase there is no drag effect and no electric current is produced by neutron vortices. For helical vortex 
motion the superfluid velocity of neutrons is 

2TO„r 

where = Vx' is the correction to the velocity due to the helical vortex motion, 

x' = iA— ^exp[i(a;i + fc||^-0)] . (17) 

In the npe phase at r < A, where A is the screening length, the induced proton velocity around the neutron 
vortex is given by (Alpar et al. 1984): 

5m* 5m* h ^ , , , , ^ 

Vp«^v„=^- — <^+v;t, 18 

m* TO* 2TO„r 

where Vp=(5TO* V^/to*. The last term in (18) is time dependent and produces density oscillation of protons 
according to the continuity equation 

^ = -p,divv;. (19) 

One should note that only z component of Vp contributes to the density oscillation since divVpj^=0. The 
Debye length for electrons, Ae = y/meC^ /^Tre^ne, is much less than characteristic scales of the proton density 

variation. As a result, free electrons screen the electric field produced by the moving protons, that is average 
electron velocity is approximately equal to the oscillating part of the proton velocity. In Aen phase there is 
no drag effect and electron motion is not coupled with neutron vortices. 

We will asstune that the core-crust interface is transparent for electrons and they can freely move from 
one phase into another. The assumption is reasonable since the electron Fermi energy is much larger than 
any Coulomb barrier that might exist at the interface over the interparticle scale. At the interface the 
oscillating electron motion generates zero sound waves which propagate across the stellar crust. To estimate 
the power of the radiated sound one can use the following simple analysis. We note that kinetic energy 
of ultrarelativistic electrons at the crust bottom is much larger than their interaction energy (Fermi gas 
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becomes more ideal with increasing its density). Therefore, the radiated power of zero sound should be 
approximately the same as in the limiting case of non interacting electrons. In this limit, the electrons move 
in a ballistic regime and the energy flux across the crust-core interface is given by 



Q = j EeVpdS, (20) 

where vp ^ cis the electron Fermi velocity and Ee is the contribution to the electron energy density due to 
the helical vortex motion. The integral in Eq. (20) is the surface integral over the interface. Eq. (20) takes 

into account that only electrons that move along the vortex line are coupled with the helical vortex motion 
and carry its energy. The rest electrons can not be excited by the helical wave because it would mean energy 
radiation by the helix in a perpendicular plane which is not possible (see Sec. 2.1). 

As we have shown above, see Eq. (10), for a vortex in a single component superfluid the change in 
kinetic energy of neutrons determines the energy of the helical wave. However, this is not the case for a 
neutron vortex in the stellar core. Due to the drag effect the helical vortex motion inevitably produces 
density oscillation of superconducting protons and normal electrons (as we have mentioned in Sec. 2.2 the 
drag is very strong). As a result, the density oscillation of electrons (which is accompanied by large change 
in their Fermi energy) gives the main contribution to the helix energy in the stellar core and, hence, the 
energy of the helical wave is equal to / E^LdS, where L is the vortex length. 

Taking into account Eq. (20), we obtain the following average energy flux of zero sound waves produced 
by a helical wave of a single vortex: 

g^-z/^-^-^^4 [exp(W)-l] - (^1) 

One should note that the wave length of a sound is much larger than the size A of the area from which it is 
generated. Therefore, the sound wave has approximately a spherical front (which is similar to diffraction of 
light on a small orifice). 

The number of vortex modes within the interval dk\\ is equal to Ldk\\/2'n = L\/2mndijj j^-K^Tiw ln(a;c/w) 
and, therefore, the spectral power density of sound waves produced by a single vortex line is 



Pr(c.)^ , , ^^^^ . (22) 

\/27rVln(a;c/w) [exp(na;/fcsT) - 1] 

Eq. (22) is a universal expression that actually depends only on the core temperature T and the neutron 

mass ?n„. Multiplying Eq. (22) by the number of neutron vortices in the superfluid A'^ = 2m„i?^r2/7i and 
dividing by 27r, we obtain the spectral density of sound waves power radiated from the oscillating vortex 
lattice in unit solid angle 

Pv(c.) = cmTR^^^ 

v^7rV^ln(a;c/a;) [exp(;ia;/A;BT) - 1] 

where T is the temperature of the NS core. 

Fast magnetosonic waves generated by vortices at the core-crust interface propagate across the stellar 
crust toward the surface. Along the wave propagation the properties of medium do not change very much 
over the length of the wave and the approximation of geometric acoustics is valid. The energy propagates 
along the rays and in the absence of absorption the energy flux remains constant along the direction of the 
ray, while properties of the medium slowly vary in space. 
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3. Attenuation of fast magnetosonic waves in the stellar crust 

In this section we estimate attenuation of fast magnetosonic waves propagating across the stellar crust. 
We assume that crust matter is everywhere in a crystal phase: T < T^, where is the melting temperature. 
For iron crust Tm — 2.5 x lO^pg^^K, where pe is the matter density in units lO^g/cm'^ (Hernquist 1984). 
Near the stellar surface there is a phase transition between a gaseous atmosphere and a condensed metallic 
phase. For iron the critical temperature of the phase transition can be estimated as Tcrit 
where H12 = which is larger than the surface temperature (~ lO^K) for strong enough magnetic 

field. We also assume that the atmosphere above the condensed surface has negligible optical depth which is 
the case when T ^ Tdn/S (Lai 2001). Under such conditions electromagnetic radiation propagates through 
the atmosphere without damping. In Appendix C we study in detail attenuation of electromagnetic waves in 
hydrogen atmosphere and show that atmosphere is transparent for electromagnetic radiation if the density 
at the atmosphere bottom p <^ 4g/cm^. The density of the metallic phase at the stellar surface depends on 
magnetic field. E.g., for iron crust and H = lO^^Gs the matter density at the crust surface is p « lO'^g/cm^ 
(Thorolfsson et al. 1998). 

In our problem the speed of fast magnetosonic wave uis larger than the electron Fermi velocity (see 
Appendix B) and the speed of usual sound in the crystal lattice (phonons) . Under such conditions the simple 
absorption processes (Landau damping) of a fast magnetosonic wave by electrons and phonons are not allowed 
due to a requirement of energy and momentum conservation. In Appendix A we discuss attenuation of zero 
sound and demonstrate that it is exponentially small for a weakly interacting Fermi gas, such as a gas 
of electrons in the NS crust and core (see Eq. (A9)). This contribution into attenuation comes from the 
collision integral in the kinetic equation. In the presence of magnetic field the electron sound is coupled 
with electromagnetic waves which results in another mechanism of damping. The point is that the fast 
magnetosonic waves are accompanied by oscillation of the electric field and in the presence of finite electrical 
conductivity a the magnetosonic waves slightly dissipate energy. This contribution dominates and we discuss 
it in this section. To find the wave attenuation one should solve equations of magnetic hydrodynamics for 
electron motion with the resistive term (see Eq. (B12)), they are derived in Appendix B. As a result we 
obtain the following dispersion equation for the magnetosonic branch 

where Us ~ vf is the speed of zero sound, ua = H/y^Anpe is the Alfven velocity, pe is the electron density 
and fc|| is the component of the wavevector along the magnetic field H. 

In the metallic crust the scattering of electrons on phonons is the dominant electron scattering mecha- 
nism. We will not discuss scattering on impurities since strong gravitational field produces space separation 
of nuclei with different masses during the "hot" stage of crust evolution. As a result, concentration of nuclei 
with different charges immersed in lattice sites is expected to be negligible. Effect of magnetic field on electron 
scattering is negligible unless the system as a whole occupies only a small number of Landau levels. Quan- 
titatively this condition means that the number of levels populated is less than about 10 (Hernquist 1984). 
For iron lattice, the condition is equivalent to 

P6 < O.eifff . (25) 

When condition (25) is satisfied the electron motion becomes effectively one dimensional. However, in ID the 
simple electron-phonon scattering is not allowed because it is impossible to satisfy equations of energy and 
momentum conservation. Therefore, at densities given by Eq. (25) the electron scattering and magnetosonic 
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attenuation is suppressed by magnetic field. The wave attenuation can be substantial at larger densities 
where many Landau levels are populated and the magnetic field does not suppress the electron scattering. 
Let us now focus on this region. In this region the electron Fermi momentum is given by pp « 1i{2>iT^n(,)^/^ 
and the Alfven velocity ua is smaller than Ug « vp (which is the case when pe > 0.16Hi2^/Z, where Z is 
the nuclear charge, A is the atomic mass number). 

At Ua <C Us we obtain from Eq. (24) the following correction to the frequency of fast magnetosonic 
wave 

i(c^ - u^)^a;^u^ sin^ 9 

- = sLy ' ^''^ 

where 6 is the angle between the wavevcctor k and the magnetic field H, Ug ^ vp is the speed of the wave. 
The wave propagates without damping a distance I which depends on electron conductivity: 

Imw ijj'^{c'^ - u1Yu\ sin^ 9 ' 
In the presence of magnetic field the electron conductivity can be estimated as (Kittel 1968) 



Re— = Re— ( 1 + ^f^"" „ ^ , Re— = Re— , (28) 

cr_L (To V +KsJ "^11 "^0 



where ctq ~ e^(?ne/ Epv^s is the static conductivity without magnetic field, u)He = eHc/Ep is the electron 
cyclotron frequency, Ep = \J ra\c^ + c^p|, is the electron Fermi energy, r^eff is the effective frequency of 
electron collisions, c7_l, cry are conductivities across and along the field H. In magnetosonic waves the 
oscillating electric field is perpendicular to the magnetic field H. Therefore, cr_L should be used to estimate 
the attenuation. 

Also in the considering density region the crust temperature is much less then the Debye temperature 
(Kittel 1968) 

Tb = = ^-"^ lO^p^/' -^K, (29) 

where Cs is the speed of usual sound, ri, is the ion concentration. In the limit T ^ Tu there are two effective 
frequencies of clcctron-phonon collisions (Lifshitz & Pitacvskii 1979). One of them i^ic,ph describes energy 
relaxation (fast relaxation) and responsible for thermal conductivity in the absence of magnetic field 



^^le,ph~^U- = „_.4/3 (30) 



The other one f2e,ph describes slow relaxation of electrons in momentum directions and determines electrical 
conductivity 

"^-'>^2^iTbj = ^ • (31) 

Using Eqs. (27), (28), (31) in the limit I'cS <C t<J <C ujhc we obtain for ultrarelativistic electrons (that is 
when pe > A/1.03Z) the following expression for the free path of fast magnetosonic waves 

^ STTc^nee^ ^ 1.14 X 103p^^/^Zi^/3 

^ (c2 - 4)2u^w2^Ef !/2e,ph sin^ 9 ~ Hf^TiA^^^ sin^ 9 ' ^ ' 
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The matter density at the stellar crust p as a function of the distance to the surface h can be obtained 
from the condition that hydrostatic pressure / pgdh is equal to the pressure of electrons (here g is the 
acceleration of gravity). The pressure of ultrarclativistic electron gas is given by 

,.i^n^r.^-^,"'iz/Ar^ (33) 

and we estimate 



Fast magnetosonic waves propagate across the stellar crust without attenuation if at any density / > h. The 
attenuation increases with decreasing the matter density. Therefore, the maximum damping occurs in the 
boundary region where magnetic field still does not suppress the electron scattering. According to Eq. (25), 
in the region of maximum damping pa « 0.6i?j2 (the region is located several meters below the surface). 
Then using Eqs. (32), (34) we obtain the following condition at which waves pass the critical region and 
reach the stellar surface without attenuation (for IgFe lattice): 

For T = lO^K, H = 4 X lO^^Gs and g = lO^^cm/s^ only waves that in the critical region propagate along 
magnetic lines within the angle |^| < 0.064rad = 3.7° (that is only about 0.2% of the total vortex radiation) 
reach the stellar surface without attenuation. However, with decrease of the surface temperature the angle 
becomes larger, e.g., if T = 5 x lO^K, we obtain \6\ < 21°. Due to anisotropic attenuation the observed vortex 
radiation should be pulsed if magnetic poles are displaced from the axis of star rotation (like in pulsars). 
The radiation is maximum when the line of sight is close to the magnetic axis. According to Eq. (35), for 
colder NSs the shape of vortex pulse becomes broader. 



4. Treinsformation of fast magnetosonic waves into electromagnetic radiation at the stellar 

surface 

Let us consider the boundary between the stellar crust and atmosphere. One can approximately treat 

the atmosphere as a vacuum because its density is much less then density of the crust. When a magnetosonic 
wave incidents on the boundary, it undergoes partial reflection and partially transforms into electromagnetic 
wave which propagates into the vacuum. The relation between the waves is determined by the boundary 
conditions at the surface of separation, which require the tangential component of the electric field and 
the normal component of the magnetic field to be equal. Conservation of the energy flux imposes another 
boundary condition (Landau & Lifshits 1960): 

= 0, (36) 

where P = Po + u^p' is the pressure, p' is the density oscillation of electron quasiparticles, and i?„ are 
tangential and normal components of the magnetic field. The square brackets denote the difference between 
the values on the two sides of the surface. From Eq. (36) we obtain the following condition for the oscillating 
part of the magnetic field h and the density p': 

= 0, [hn] = 0. (37) 



P + 



uIp' + 
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The problem can be substantially simplified because near the stellar surface the Alfven velocity is 
much larger then the speed of sound ua ^ Us- For H = lO^^Gs, pe — 10^{me/mp)g/cm'^ we obtain 
itA = 4 X lO-'^^cm/s ^ c. In this case the second term in Eq. (37) is in the factor u\/ul ^ 1 greater than the 
pressure oscillation u^p' and the boundary condition reduces to [hf] = 0. Also, under the condition ua ^ Ug 
the fast magnetosonic wave propagates with the velocity u = ua/ + uj^/c^ which is independent of the 
wavevector k (see Eq. (B23)). The particle velocity V in the wave lies in the kH plane and perpendicular 
to H, while oscillating magnetic h and electric fields E are perpendicular to each other and k: 

E= ^^ {hxk), (38) 

where fc is a unit vector in the direction of k (see Fig. 3). We obtain that the fast magnetosonic wave 
near the stellar surface is similar to an electromagnetic wave in a medium with the dielectric constant 
£ = 1 + /u\ and the boundary conditions at the interface coincide with those for electric and magnetic 
field at the boundary of the dielectric medium. Therefore, transformation of the fast magnetosonic wave into 
electromagnetic one is described by the same formulas as refraction of electromagnetic waves at the interface 
between the dielectric and vacuum. 

For example, let us consider the case when the magnetic field H lies in the plane of incident wave. Then, 
the oscillating electric field E is perpendicular to the plane of incidence. The transmission coefiicient is given 

by 

sin(2(9i)sin(26'2) 
sin2(^i+^2) 

where 6i is the angle of incidence and 62 is the angle of refraction (sin 62 = \/e sin ^1 ) . For normal incidence 
(^1 = 0) the transmission coefficient is 

4Vi _ 4v/l + cVu^ ^^^^ 



If » c (which is the case when H > lO^^Gs) the magnetosonic wave is completely transformed into 
electromagnetic wave (r w 1). In the opposite limit ua ^ c the transmission coefiicient is small and r 
K, Aua/c. Fig. 4 shows the general scheme of the vortex mechanism of NS radiation. 



5. Spectrum of neutron star radiation 

The spectral density of sound waves power radiated from the stellar core in unit solid angle is given 
by Eq. (23). However, as we have shown, only a fraction of magnetosonic waves that in the critical region 

propagate in the direction close to the magnetic lines reaches the stellar surface and only a fraction is 
transformed into electromagnetic waves which propagate into surrounding space. The fraction depends, in 
particular, on magnetic field strength and its distribution over the stellar surface. 

Now let us discuss the region of applicability of Eq. (23). There is a fundamental limitation that 
affects the spectrum of vortex radiation in the infrared band. The limitation is that at any frequency the 
electromagnetic radiation produced by vortices should not exceed black body radiation produced by the 
surface AttB? with the same temperature T. This statement is one of the formulations of Kirchhoff's law 
which is a consequence of general thermodynamic principles. We assume that vortices are excited thermally 
and if at some frequency the vortex radiation becomes comparable with radiation of a black body this means 
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that at such frequency the rate of thermal excitation imposes the main restriction on the radiated power: 
vortices radiate maximum power which can be pumped from the thermal reservoir. Kirchhoff 's law results 
in the following limitation on the frequency at which Eq. (23) describes the spectrum of electromagnetic 
radiation produced by vortices: 

27r > ?i3/5 1^1/5 ' ( ) 

where lOc = ?i/4m„^^ « 10^°Hz. For Q,/2'!t = 4s~^ we obtain a;/27r > 1.5 x lO^'^Rz. At lower frequencies the 
spectrum of vortex radiation (in unit solid angle) follows Planck's formula 

~ 47r2c2 [exp(;ja;/fcBT) - 1] ~ Att^c' ' ^ ' 

where T is the temperature of the NS core. 

One should mention that Eq. (41) for the kink frequency contains no free parameters. To take into 
account the gravitation redshift z one should divide Eq. (41) by 1 + z. Consequently, the observation of the 
kink in the infrared band allows measurement of the gravitation redshift at the core surface. 

To compare our theory with available spectrum observations in different frequency bands it is convenient 
to represent the spectral density of vortex radiation Pv(w) as a power law Pv(w) oc w", where a is the spectral 
index. The ^J\n{uJ^uJ) function in the denominator of Eq. (23) shifts the spectral index by a small value 
l/2hi{uJc/Lij). The shift depends weakly on lo and changes the spectral index of vortex radiation from 
a « —0.46 in the optical band to a w —0.43 in the X— ray band. 

Apart from vortex contribution there is thermal radiation from the NS surface, which for middle-aged 
pulsars dominates in the ultraviolet and soft X— ray bands. One should note that spectrum of thermal 
radiation may substantially differ from a black body. E. g., this is the ease for NS atmospheres composed 
from light elements (hydrogen, helium) (Romani 1987). However, despite of the big difference, the observed 
thermal spectrum can be well fitted with a hydrogen (helium) atmosphere model as well as with a black body 
spectrum (Yakovlev et al. 1999a). The point is that spectrum of thermal radiation is usually detected only in 
the soft X— ray band. No data are available in the ultraviolet band due to strong interstellar absorption. In 
the soft X— ray band one can successfully fit the data assuming two black body components: hot polar caps 
and contribution from the rest (colder) NS surface. In the case of hydrogen (helium) atmosphere the data 
are well fitted with only one component contribution. In the optical-ultraviolet band the two models predict 
substantially different spectra and the hydrogen (helium) atmosphere curve lies above the black body giving 
the intensity several times greater. However due to the lack of data in this spectral region it is difficult to 
choose the adequate atmosphere model. 

In this paper we estimate thermal radiation assuming helium atmosphere and use the numerical results 
obtained by Romani (Romani 1987). The model fits well the thermal spectrum with only one free parameter 
(the effective surface temperature Tco). However, the main reason for choice of helium atmosphere is that if 
there is vortex radiation from the stellar surface then the atmosphere, in principle, can not be a black body 
because black body surface layer would absorb all passing radiation. From this point of view the detection 
of vortex radiation is a strong argument against black body atmosphere models. 

The total radiation from a NS is given by the sum of vortex (Eqs. (23), (42)) and thermal components: 

P(a;) = aPv(w) + -Pth(w), (43) 

where Pth(<^) is the spectrum of surface thermal radiation. In Eq. (43) we introduced a dimensionless 
parameter a. The parameter is useful when we compare theoretical predictions with radiation observed from 
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pulsars. The free parameter a takes into account partial absorption of the vortex radiation near the stellar 
surface, geometrical effects related to unknown magnetic field distribution and position of the line of site 
with respect to the cone swept by star's magnetic axis. According to Eq. (35), for typical parameters of NSs 
and the surface temperature Tg rs lO^K only about 0.2% of vortex radiation reaches the stellar surface and 
transforms into electromagnetic waves. However, the outgoing radiation is concentrated in a small solid angle 
near the magnetic axis which is determined by distribution of magnetic field. For typical pulse shapes this 
angle is of the order of O.lrad. That is parameter a should be of the order of 0.002 • 27r/0.1 w 0.1. However, 
according to Eq. (35), with decreasing the surface temperature the parameter a increases as a oc l/T^ and 
for colder NSs becomes of the order of 1. 

In next section wc compare the vortex radiation for different NSs with the thermal surface contribution. 
In such relative comparison the neutron star radius Rg (which is approximately equal to the core radius R) 
cancels in equations. As a result, in estimates there are only two free parameters: the core temperature T 
and the geometrical factor a. The effective surface temperature T^n is directly determined from the thermal 
spectrum. 

6. Discussion 

Since the discovery of pulsars in 1967 numerous attempts were made to detect radiation from a NS 
surface. Observations were mainly successful in soft X— ray band because of high thermal luminosity in this 
region. Recent observations with ROSAT (0.24 - 5.8 x lO^Hz), ASCA (0.17 - 2.4 x lO^^Hz) and Chandra 
obtained a detailed spectral information of a number of pulsars at X-ray energies. Only a few pulsars were 
observed in optical and ultraviolet bands and only some of them have detectable emission in these bands. In 
Fig. 5 we compare the observed radiation spectrum of middle-aged pulsars PSR B0656-I-14 {Q/2n = 2.6s~^) 
and Vela (17/27r = 11.2s~"'^) with the radiation spectrum predicted by our theory. Thermal radiation of the 
stellar surface dominates in the ultraviolet and soft X— ray bands. We assume helium atmosphere with the 
effective surface temperature T^s = 4.9 x lO^K for PSR B0656+14 and T^s = 7.8 x lO^K for Vela. The 
NS radius i?,, is related to the pulsar distance D: i?J10km = 5.28\/l + z{D/lkpc) for PSR B0656+14 and 
Rs/lOkm = 4.lVTTI(£'/lkpc) for Vela (z is the NS rcdshift). For PSR B0656+14 and D = 0.25kpc we 
obtain Rg « + zkm, while for Vela and D = 0.35kpc the surface radius is Rg ~ 14V1 + zkm. The 

vortex contribution (dash line) dominates in the infrared, optical and hard X— ray bands, where its spectrum 
has a slope a « —0.45. In the far infrared band vortices radiate the maximum power which can be pumped 
by thermal excitation. In this band the radiation spectrum changes its behavior and follows Planck's formula 
with P{lo) oc u)^. The sum of the vortex and thermal components is displayed with the solid line in Fig. 
5. For PSR B0656-I-14 we take the core temperature T = 6.4 x lO^K and the geometrical factor a = 0.18, 
while for Vela T = 8 x lO^K and a = 0.33. The observed broad-band spectrum is consistent with our 
model for typical NS parameters. This suggests that the vortex mechanism of radiation operates in a broad 
frequency range from IR to hard X— rays. One should mention that the measured IR-optical spectrum has 
nonmonotonic behavior, which probably indicates the presence of unresolved spectral features or could be 
partially due to contamination of the pulsar flux by the nearby extended objects (Koptsevich et al. 2001). 
The slope of the non-thermal component is less steep than it was deduced by Pavlov et al. 1997 and Kurt et 
al. 1998 based on the optical-UV data only. Also, in the optical band the radiation of middle-aged pulsars 
was found to be highly pulsed with pulses similar to those in the hard X— ray band (Shearer et al. 1998; 
Shearer et al. 1997). This agrees with the vortex mechanism which predicts pulsed radiation. 

In Fig. 6 we plot the broadband spectrum of the Geminga pulsar (0/27r = 4.2s~^) in optical— X- 
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ray bands. Solid line is the fit by the sum of the vortex and thermal contributions. In calculations we 
take Tcff = 2 x IQ-'^K, T = 3 x lO^K and a ^ 0.22. The radius of the Geminga pulsar is: i^^/lOkm = 
5.21vT+T(D/lkpc), which for D = 0.16kpc gives Rs ~ 8.3V1 + ^km. Our fit reasonably agrees with the 
observed spectrum apart from the dips in the IR-optical bands which probably appear due to ion cyclotron 
absorption (Martin ct al. 1998). Also there is strong scattering of the data in the hard X— ray band which 
makes difficult comparison with the theory and, hence, estimates are less reliable. The data scattering might 
indicate on unresolved emission and absorption lines in this spectral region. 

It is worth to note that the formula for vortex radiation contains two free parameters (T and a). As 
a result we can not determine each of these parameters separately from the observed fluxes, only their 
combination. However, the spectral index of vortex radiation is a fixed parameter in cnir theory. If a single 
mechanism of nonthermal radiation operates in the broad range, then the continuation of the X— ray fit to 
the optical range determines the spectral index of the nonthermal component with a big accuracy. Good 
quantitative agreement of our theory with the observed spectral index serves as a strong evidence that 
the vortex mechanism is responsible for radiation of middle-aged NSs in the IR, optical and hard X— ray 
bands. Also, there is an evidence of vortex contribution in radiation of young Crab pulsar. In optical band 
an emission during the "bridge" phase interval has been detected from the Crab with the spectral index 
dopt = —0.44 ± 0.19 (Golden et al. 2000). The index differs from the flat spectrum of the peaks and can 
indicate on the vortex component. 

Moreover, according to Eq. (23), the vortex contribution exponentially decreases at Tilo > ksT. Obser- 
vation of such spectrum behavior allows us to determine directly the temperature of the stellar core T and 
can be a possible test of our mechanism. For PSR B0656+14 and the Vela pulsars the power law spectrum 
in the hard X— ray band shows no changes up to the highest frequency 2 x lO^^Hz at which the data are 
available (see Fig. 5). According to our theory, this indicates that temperature of the NSs core is larger than 
2 X 10*K. Measurements of spectra of middle-aged pulsars in the 10^^-lO^^Hz range are needed to search for 
possible manifestations of the core temperature. 

One should note that recent observation of the Vela pulsar with the Rossi X— ray Timing Explorer has 
covered the energy band 2 — 30keV (0.49 — 7.3 x lO-'^^Hz). These data in combination with OSSE observations 
(70 — 570kcV) allows us to estimate the temperature of the Vela core. Light curves of the Vela pulsar have 
several peaks. We associate the vortex radiation with the second optical peak (see Fig. 1 in Harding et 
al. (1999) and Fig. 4 in Strickman et al. (1999)). In the ROSAT band (0.06 - 2.4keV) the vortex peak is 
weakly distinguishable on the background of strong thermal emission. However, the peak is clearly seen in 
the Rossi band in which thermal radiation is negligible. The peak becomes very small in the OSSE band 
and completely disappears at EGRET energies (> lOOMeV). The vortex peak gives the main contribution 
to the average radiation intensity in the optical and hard X-ray bands and its spectral index agrees with 
our theory. One should note that in the Rossi band other peaks appear in the Vela light curves (Harding 
et al. 2002). Position of those peaks coincides with pulses at EGRET energies and, therefore, they have a 
different origin. The new peaks give the main contribution to the average radiation intensity at energy greater 
then lOkeV. As a result, at such energies the average intensity carries no information about vortex radiation 
and a light curves study is necessary to separate the vortex contribTition. To estimate the core temperature 
one should trace the evolution of the vortex peak in the energy interval IkeV— IMeV. According to Eq. (23), 
the frequency dependence of the peak amplitude depends only on the core temperature T. To our knowledge 
such analysis of the vortex peak amplitude has not been done yet. Nevertheless, it is possible to estimate 
the core temperature from available data. If we assume that the amplitude of the vortex peak follows Eq. 
(23), then we can find total (integrated over frequency) radiation intensities of vortices in different energy 
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bands. Then wc compare the detected integrated intensities of the vortex peak in the Rossi and OSSE bands 
(Strickman et al. 1999) with those predicted by Eq. (23). Wc found that if there was no temperature decay 
of the vortex peak the total vortex intensity in the OSSE band should be three times larger than those 
actually observed. This indicates on the decay of the vortex spectrum and results in the following estimate 
of the core temperature of the Vela pulsar: T w 8 x lO^K. 

The temperature estimate (T > lO^K) allows us to make a conclusion about the interior constitution of 
the NSs. A hot NS cools mainly via neutrino emission from its core. Neutrino emission rates, and hence the 
temperature of the central part of a NS, depend on the properties of dense matter. If the direct Urea process 
for nucleons {n ^ p + e + Ug, P + e ^ n + Ve) is allowed, the characteristic cooling time is Turca ~ 1 min /Tg, 
where Tg is the final temperature measured in imits of lO^K (Pethick 1992). Such neutron star cools to 
lO^K in weeks. However, the characteristic age of PSR B0656+14 is about lO^yrs, while the Vela pulsar is 
about 10^ years old. So, the NS core cools down to lO^K at least 10^ times slower then the rate predicted 
by the direct Urea process. This estimate excludes equation of states incorporating Bose condensations of 
pions or kaons and quark matter. The point is that the neutrino emission processes for these states may 
be regarded as variants of the direct Urea process for nucleons (Pethick 1992). As a result, all these states 
give rise to neutrino emission, though generally smaller, comparable to that from the direct Urea process 
for nucleons which is inconsistent with the discrepancy of the cooling rate in the factor 10^. The conclusion 
about absence of exotic states of matter in the stellar core agrees with recent measurement of gravitational 
redshift from a NS (Cottam et al. 2002). 

It is worth to mention here that dissipation of rotation kinetic energy in principle can produce heating 
of the stellar interior and supply energy for neutrino emission (Riper et al. 1995; Larson & Link 1999). Such 
possibility correlates with recent observation of emission properties of 27 pulsars by Becker and Triimper 
(1997) which suggests that pulsars radiation is emitted at the expense of their rotational energy (for ^1/2^ = 
10s~^ the rotation energy ~ MqR'^Q^ fa 10^"erg). The superfluid core possesses the main part of the 
rotation kinetic energy, while radiation and dissipation of the angular momentum occurs from the stellar 
surface. As a result, the outer crust slows down faster than the superfluid core. The difference in the angular 
velocities between the superfluid and normal components causes friction which heats the core and transfers 
the angular momentum from the interior part to the outer crust. Superfluid vortices can play a crucial role in 
such process (Zhuravlev et al. 2001). One should note, however, that the neutrino emission rate comparable 
to that from the direct Urea process also contradicts to the rotation energy loss of NSs. Indeed, according 
to Lattimer et al. (1991), the neutrino emissivity in the direct Urea process is Qurca ~ 10~^^T^erg/s • cm^ 
and the total power radiated by neutrino at T = 10*K is Purca = l^ri^^Qurca ~ 4 x lO^^erg/s. This value is 
much larger then the rate of the kinetic energy loss of the middle-aged pulsars 10^^ — lO^^erg/s. 

One should note that the characteristic time for cooling by the modified Urea process for normal nucleons 

is TmodUrca ^ ^Y'c/Tq (Pcthick 1992). That is if the direct Urea process cannot occur, the core temperature 
will exceed ~ lO^K for ~ 10^ years. This result does not contradict the estimate T > lO^K for middle-aged 
pulsars. However, more accurate determination of the core temperature of the Vela pulsar (T w 8 x lO^K) 
demonstrates that NSs cool down much more slowly than could be explained by the modified Urea process for 
normal nucleons (such process would cool the core down to 8 x lO^K during 4yrs, however the age of the Vela 
pulsar is about lO^yrs). Possible explanation of this observation is that superfluidities of neutrons and protons 
in the stellar core are very strong so they substantially reduce all neutrino processes involving nucleons. If 
ksT <C A, where A is the gap in the excitation spectrum of a superfluid, the modified Urea rates are reduced 
by a factor ^ exp(— 2A/A:bT). Using the measured value of the core temperature one can estimate the critical 
temperature of the superfluid transition. Theory predicts that in the stellar core the critical temperature of 
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the superconducting transition for protons Tcp is larger than the transition temperature of triplet neutron 
superfluidity Ten (Hoffberg et al. 1970; Amundsen & Ostgaard 1985a; Amundsen & Ostgaard 1985b). As a 
result, the proton gap gives the main contribution to the suppression of neutrino emission by the modified 
Urea process. The corresponding reduction factors can be found in (Yakovlev et al. 2001). However, the 
appearance of neutron or proton superfluid initiates an additional specific neutrino production mechanism 
due to Cooper pair formation (Yakovlev et al. 1999b). Such mechanism of neutrino emission dominates until it 
dies out at T < O.lTc. The process involving superfluid neutrons n — > n+v+D is more important in the stellar 
core since their number greatly exceeds the number of protons. To find Ten we assume triplet-state pairing 
of neutrons with the gap A = Ao\/l + 3 cos^ 9. The neutrino emissivity due to Cooper pairing of neutrons 
can be estimated using the result of Yakovlev et al. (1999b) Qcp = 2.2 x 10-'^22^^i^(r/rc„)erg/s • cm^, 
where F{x) w 1.27exp(— 2.376/x)/a;^ for a; <C 1. The stellar core cools down to the temperature T during 
the time r = dTC{T)/Qcp{T), where C(T) is the specific heat. When T < T,, the neutron (proton) 
contribution to the specific heat is exponentially small. Therefore, C(T) is determined by the electron 
component C(T) = {^ir^Y/^n^/^k'^T /2ich, where n « lO'^^cm"'^ is the electron concentration in the stellar 
core. Using parameters of the Vela pulsar T w 8x lO^K, r w lO^yrs, we find the following critical temperature 
of triplet neutron superfluidity in the center part of the Vela core Ten ^ 9 x lO^K. Again one should mention 
the possibility of heating the stellar core by dissipation of the rotation energy. If such process is responsible 
for the Vela temperature than in equilibrium the heating rate is equal to the total power radiated by neutrino 
(4/3)7ri?^(5cp- The kinetic energy loss of the Vela pulsar is 6.9 x lO^^erg/s (Becker & Triimper 1997). If we 
assume that all kinetic energy loss goes in the neutrino emission than we obtain Ten ~ 6 x lO^K {R k, IQkm). 
The true value of Ten fies between this estimate and the previous value 9x lO^K found under the assumption of 
no extra heating. More detailed analysis shows that the core temperature of middle-aged pulsars is probably 
maintained by the heating due to the dissipation of rotation energy. Our estimate of the critical temperature 
of neutron superfluidity agrees well with realistic nucleon-nucleon interaction potentials that fit recent world 
data for pp and np scattering up to 500Mev (Baldo et al. 1998). Such potentials give Ten ~ 7.5 x lO^K 
at the density 2'^o that is expected at the center region of a canonical neutron star without exotic inner 
core, and predict Ten ~ 2 x lO^K at the core-crust interface. One should mention that precise knowledge 
about nucleon-nucleon interaction at high energies (up to ^ 2.5Gcv) is a result of progress in accelerator 
physics during the last decade (Arndt et al. 1994; Arndt et al. 1997). Estimates of Ten made in 1970's 
and 1980's usually substantially differ from the current value. E.g., Amundsen & Ostgaard (1985b) gave 
the maximum value of Ten in the stellar core Ten = 9 x lO^K which is almost one ordc^r of magnitude less 
then the current result. Recent experiments on nucleon-nucleon scattering favor in strong superfiuidity and 
superconductivity in the stellar core, such observations are also supported by the theory (Elgaroy et al. 1996; 
Baldo & Grasso 2001). 

It is worth to note that we may rule out the exotic states of matter only for NSs which display vortex 
radiation in the hard ray band. However it may not be the case for the whole NS population and a 
systematic study of pulsars radiation in the X— ray band is necessary. If exotic states do exist in some dense 
NSs the middle-aged pulsars with the core temperature < lO^K could be candidates for stars with Bose 
condensates or quark matter. 

7. Energy transport across the stellar crust 

The theory of NS cooling has been developing over more than 30 years. First let us discuss the current 
understanding of temperature distribution in a crust of middle-aged NSs. The question is important because 
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NS cooling theory is used as a method to determine the core temperature from the surface temperature 
and, hence, to explore properties of superdense matter. In current theories it is assumed that high density 
of the stellar interior provides high thermal conductivity of internal layers. As a result, in about 10^ — 10^ 
years after the star formation, a wide, almost isothermal region is formed within the entire core and the 
main fraction of the stellar crust (Yakovlev et al. 1999a). The isothermal stellar interior is surrounded by a 
very thin subphotospheric layer with thickness of the order of several meters. This low dense surface layer 
produces thermal isolation of the NS interior. The heat transport through the isolating layer determines the 
relation between the surface temperature Tg and the interior temperature T. In particular, T/Ts depends 
on the surface gravity and the chemical composition of the outer NS envelope. For Ts = 5 x lO^K and a 
nonmagnctizcd NS the ratio T/Tg typically lies in the interval 20 100 (Potckhin ct al. 1997). 

Possible presence of a strong magnetic field increases the uncertainty. The point is that magnetic field 
suppresses the electron thermal conductivity across the field lines and enhances along the lines. The effect of 
magnetic field on NS cooling has been studied by several authors. E.g., a relationship between the internal 

and surface temperatures in a NS with the magnetic field normal to the surface was investigated by Van 
Riper (1988). Page (1995) as well as Shibanov & Yakovlev (1996) considered the NS cooling with a dipole 
magnetic field. 

As we have found based on spectra observations the effective surface temperature of the Vela Pulsar is 
7.8 X lO^K, while the core temperature is about 8 x lO^K. This value of the core temperature is much larger 

than those predicted in heat transport models discussed in literature so far. The purpose of this section is 
to explain the big difference. Here we show that if near the core surface the magnetic field is parallel to the 
core it produces substantial heat isolation and can explain the big difference between the core and surface 
temperatures. Such assumption is simply a boundary condition for the magnetic field at the superconductor 
boundary. In this model the main temperature drop occurs in a thin layer near the crust bottom which 
differs from current cooling theories that predict isothermal inner crust. 

Let us consider heat transport across the stellar crust in the presence of a strong magnetic field. Wc 
assume that free electrons give the main contribution to thermal conductivity. In the presence of a strong 
magnetic field {Ep '> WHe ^ t^e) the electron thermal conductivity across the field can be estimated as 

where is the frequency of electron relaxation in momentum directions which can be due to electron-phonon 
(Eq. (31)) or electron-electron collisions. At T ^ Tp, where Tp is the Fermi temperature, the frequency of 
electron-electron collisions for ultrarelativistic electrons is 

For rie = lO'^^cm^'', T = 10*K we obtain v^.c ~ 3.2 x 10^°s^^. Near the crust-core interface i/c,c > J^2o.ph 
and, therefore, should be substituted into Eq. (44). As a result, the electron thermal conductivity across 
the magnetic field is given by 

ACex « 8.5 X 10-29 "ejl erg ^^g^ 



For Ue = lO^^cm-^, T = lO^K, H = W^Gs we obtain « 8.5 x lO^^erg/K • cm • s. Magnetic field 
suppresses /te_L in a factor i'e,e'^ie,ph/'^He ~ 10~^^. One should mention that Eq. (46) is valid when ^ 
i^e,e*^ie,ph- In the Opposite limit the magnetic field does not affect the thermal conduction. 
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Eq. (46) has to be compared with the phonon contribution to thermal conductivity which for T <C To 
can be estimated as (for phonon relaxation on electrons): 

5?i c K • cm • s 

For T = lO^K we obtain Kph « 2 x lO^^erg/K • cm • s. At such temperatures, the electron thermoconductivity 
Ke_L is still larger than Kph- 

Now let us show that magnetic thermal isolation can substantially reduce the heat transport from the 
stellar core. Under such condition the magnetosonic waves produced by vortices can become the dominant 
mechanism of energy transfer to the stellar surface. The superconducting core expels magnetic field from 
its interior so that outside the core the magnetic lines are parallel to the superconductor boundary. As a 
result, near the crust bottom the magnetic field is tangential to the stellar core which produces substantial 
thermal isolation. To estimate the heat transfer across the magnetic layer we consider the following model. 
We assume that there is a strong magnetic field parallel to the core surface in a thin layer of width d. The 
width d is much less then the crust thickness {d ^ 1km). Outside the layer we assume the magnetic field 
to be weak or not tangential to the surface so that outside the layer the field does not suppress thermal 
conductivity across the crust. 

One should also take into account that at high temperature the neutrino emission processes result in 
substantial energy loss and can affect the temperature profile in the inner crust. We assume that neutrino 
emission due to Cooper pairing of neutrons is the main mechanism of neutrino generation near the crust 
bottom. The process represents neutrino pair emission (of any flavor) by a neutron whose dispersion relation 
contains an energy gap: n — > n+u+P (Yakovlev et al. 1999b). In this process the main neutrino energy release 
takes place in the temperature range 0.2Tc <, T <^ 0.96Tc, where Tc is the critical temperature of neutron 
superfluidity. At densities near the crust bottom Tc « lO^K and, therefore, the neutrino emission due to 
Cooper pairing of neutrons modifies the temperature distribution at 2 x 10^ ^ lO^K. In this temperature 
range we assume the following neutrino emissivity Q « lO^^Tgcrg/s • cm^ (Yakovlev et al. 1999b). 

Fig. 7 shows the temperature distribution in the NS crust as a function of the distance to the crust-core 
interface z. We take temperature at the crust-core interface T = 8 x lO^K (the Vela pulsar), the width of 
magnetic layer is c? = 50m. Thermal conductivity across the layer is given by Eq. (46) with rie = 5 x lO^^cm"^ 
and H — Ax lO^^Gs. Neutrino emission substantially modifies the temperature distribution and reduces the 
heat flux at 2; < 40m. At larger distances from the core the temperature becomes low enough so that neutrino 
emission is no longer important. In this region the heat flux is conserved. At z > d the thermal conductivity 
is not suppressed by magnetic field and temperature remains approximately constant. Near the star surface 
(at distance less then Im from the surface) there is another temperature drop due to substantial decrease 
of thermal conductivity at low densities. In estimates we take the crust width is equal to 1km. Practically 
all heat flux from the stellar core is lost by neutrino emission in a thin layer near the crust bottom. In our 
model the total energy losses of the stellar core due to thermoconductivity is 7 x lO'^^erg/s (for core radius 
R = 10km), while only 1.4 x lO^^erg/s passes the magnetic layer and contribute to thermal radiation from 
the surface. 

One should note that this value is about 20 times less then thermal energy actually radiated from 
the Vela surface: for = 10km, Tos = 7.8 x lO^K we obtain AnR^aT^^^ = 2.6 x 10^2(.j.g/s g^^h big 
discrepancy allows us to advance a hypothesis that vortex radiation, rather then thermoconductivity, is 
the main mechanism of heat transfer to the stellar surface. To check the hypothesis let us estimate the 
temperature of a black body that radiates the same total power as radiated by vortices from the superfluid 
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core. Using Eq. (23) we obtain 



k]i\n^'^{huo,/kBT) ■ 



(48) 



In our model the fast magnetosonic waves generated by vortices propagate practically across the whole 
crust in a ballistic regime. For not too cold NSs a thin layer near the stellar surface absorbs practically 
all magnetosonic waves apart from a small fraction which reaches the surface. If thermoconductivity across 
the crust is suppressed by magnetic field, then vortex radiation can become the main mechanism of energy 
transfer from the core to the surface. Under such condition the energy radiated by vortices should be 
equal to the thermal energy radiated from the stellar surface and Eq. (48) determines the effective surface 
temperature which is directly measured from observations in the soft X— ray band. One should note that 
Eq. (48) contains no free parameters because T^s and T can be measured independently and an uncertainty 
in (jjc = h/AninS,'^ under the logarithm might change the answer only in a few percents. 

From observation of the Vela pulsar in the hard X— ray and soft 7— ray bands we estimated the Vela 
core temperature T 8 x lO^K. For this value of the core temperature, 0/27r = 11.2s~^ and Wc/^tt = lO^^Hz 
Eq. (48) gives T^g = 7.6 x lO^K, which approximately coincides with the effective surface temperature of 
the Vela pulsar obtained from observations of thermal spectrum in the soft X— ray band T^s — 7.8 x lO^K. 
Such excellent agreement with observations is remarkable because the theory contains no free parameters! 
This confirms our hypothesis that vortex radiation, rather than thermoconductivity, is the main mechanism 
of heat transfer to the stellar surface in middle-aged NSs (at least in some of them). Also such agreement is 
a strong argument in favor of the theory of vortex radiation. 

As a result, one can use Eq. (48) to estimate the core temperature T of middle-aged NSs from their 
surface temperature T^s and obtain the following practical relation: 



ForPSRB0656+14rcff = 4.9xlO-'^K, f^/27r = 2.6s"i and Eq. (49) gives the core temperature T « 6.4 x lO^K. 
For Geminga T^g = 2 x lO^K, Q/2Tr = 4.2s"^and we obtain T w 0.3 x lO^K. One should mention that for 
colder NSs (such as Geminga) absorption of vortex radiation near the surface becomes small (see Eq. (35)). 
In this case Eq. (49) is less reliable and can substantially underestimate the core temperature. For colder 
NSs with still strong magnetic field the total power of vortex radiation becomes larger than thermal surface 
emission. 



In this paper we discuss a new possible mechanism of NS radiation produced by vortices in the superfluid 
core. The key idea is that neutron star interior is transparent for electron zero sound, the same way as it 
is transparent for neutrinos. The waves of zero sound (in electron liquid) are generated at the core-crust 
interface by superfluid vortices which undergo thermal helical oscillations. Sound waves produced by vortices 
propagate across the stellar crust towards the surface. In the presence of magnetic field sound is coupled 
with electromagnetic waves. Virtually, they are limiting cases of the same essence - the fast magnetosonic 
wave which is known in magnetic hydrodynamics. Near the crust bottom, where the density is large, the fast 
magnetosonic waves reduce to zero sound (the effect of magnetic field is negligible). However, near the stellar 
surface the Alfven velocity becomes larger than the speed of sound and the fast magnetosonic waves behave as 



Tgln^/^ (48/T8) 



(f^/27r)2/3- 



(49) 



8. Conclusions and outlook 
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electromagnetic waves in a medium. As a result, at the stellar surfac;c the fast magnetosonic waves transform 
into electromagnetic radiation in vacuum the same way as refraction of usual electromagnetic waves at the 
dielectric- vacuum interface. For typical magnetic fields of the middle-aged NSs the transformation coefficient 
is close to 1. 

Finite conductivity of the crust results in partial absorption of magnetosonic waves in a thin layer 
near the stellar surface. As a result, only a fraction of waves which in the; absorbing layer propagate in the 
direction close to magnetic lines reach the surface and transform into electromagnetic radiation. The fraction 
depends on the magnetic field strength and surface temperature. For typical parameters of middle aged NSs 
the fraction is less than 1%, however for colder stars, such as Geminga, it can be substantially larger (even 
of the order of 100%). As a consequence of anisotropic attenuation and star rotation the observed vortex 
radiation is pulsed, pulse shape is broader for colder stars. The part of vortex radiation absorbed near the 
surface contributes to the heat transfer from the core to the surface. Estimates show that at least for some 
NSs such mechanism, rather than thermoconductivity, is the main source of heating the stellar surface. 

Stellar radiation produced by vortices has the spectral index a « —0.45 and for middle-aged pulsars 
dominates in the IR, optical and hard ray bands. While the thermal surface radiation contributes mainly 
in the ultraviolet and soft X— ray bands. The observed spectra of PSR B0656-I-14, Vela and Geminga pulsars 
agree with our theory. In particular, there is an excellent agreement with observations in two predictions 
where the theory contains no free parameters: the spectral index of vortex radiation a and the relation 
between the core and surface temperatures (given by Eq. (48)). There are also several predictions which 
have to be confirmed by future observations. For example, the theory predicts a kink of the vortex radiation 
spectrum in the infrared band (see Eqs. (41), (42) and Figs. 5,6) and an exponential decay of the spectrum 
in the hard X— ray band at Tioj > ksT (see Eq. (23)) which was not measured in details yet. 

The vortex radiation is produced in the stellar interior and contains important information about prop- 
erties of the NS core, in particular, it allows direct determination of the core temperature. Comparing theory 
with available spectra observations we found that the core temperature of the Vela pulsar is T « 8 x lO^K, 
which is the first measurement of the temperature of a neutron star core. It is important to note that the 
temperature we found is much larger than those predicted in heat transport models discussed in literature 
so far. To understand the big difference we consider heat transport in a model with magnetic field parallel 
to the core surface (which is just a boundary condition at the superconductor boundary) and show that the 
effect of magnetic thermal isolation can explain so high value of the core temperature as compared with the 
stellar surface temperature. 

A direct measurement of the core temperature by means of detection the vortex radiation opens a new 
perspective in study superdense matter. In particular, our estimate of the interior temperatures of PSR 
B0656-I-14, Vela and Geminga pulsars rules out possibility of exotic states, such as Bose condensation of 
pions or kaons and quark matter, in these objects. Also it allows us to extract the superfluid transition 
temperature for neutrons in the center region of the stellar core Tc = (7.5 ±1.5) x lO^K which is the first 
determination of such parameter from observations. 

Detection of vortex radiation opens a possibility to study composition of NS crust (NS interior spec- 
troscopy). Since magnetosonic waves generated by vortices propagate through the stellar interior the 
spectrum of vortex radiation should contain (redshifted) absorption lines which correspond to low energy 
excitations of nuclei that form NS crust. E. g., the ^^Pe nucleus has an excited state with the energy 
14.4kev= 3.5 x lO^^Hz, which would produce an absorption line in the hard X— ray band of vortex radiation 
spectrum if temperature of the stellar core is greater than 1.7 x lO^K. The other possible absorption lines 
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in the X ray band can be produced by nuclei 2iSc (12.4kev), fsMn (26kev), IfGe (13.5kev) or by heavy 
nuclei gf^U (13kev), gy^Bk (8.8kev). One should mention that bottom layers of neutron star crust may 
contain exotic nuclei with mass number up to about 600 and vortex radiation can be a possible tool to 
study their properties. However, the absorption lines of nuclei are narrow and large spectral resolution is 
necessary for their detection. If temperature of the stellar crust Tcmst ^ lO^K the relative thermal line 
width is \/2kBTcrust/Mc^ « 2 x 10~^. Under such condition the star rotation gives the main contribution 
to the line broadening and the corresponding Doppler line width for middle-aged pulsars can be estimated 
as nRs/c = 0.001 - 0.01. 

Another challenging goal for future study is interaction of zero sound with exotic states of matter which 
might exist in more dense NSs. The point is that generation of zero sound by vortices is only one of the 
possible mechanisms. If an exotic state of matter absorbs zero sound at some frequency, then, according 
to Kirchhoff's law, it will radiate zero sound at this frequency. As a result, spectrum of stellar radiation 
must contain characteristic emission lines corresponding to such processes which allows a direct spectroscopic 
study of superdense matter. In the limiting case when the inner core matter absorbs all zero sound waves 
the stellar radiation spectrum should be close to radiation of a black body with the temperature of the core. 
One should look for such objects among bright galactic X— ray sources. Also zero sound can be emitted in 
reactions between elementary particles that constitute the stellar core. It is known that some pulsars are 
powerful sources of 7— rays and it would be interesting to study the connection between 7 radiation and 
zero sound. It is also interesting to investigate a possible connection between electron sound and nonthermal 
radiation of young (Crab like) pulsars which reveals flat spectrum in the optical band. 

I am very grateful to A. Fetter, G. Shlyapnikov, M. Binger, S. T. Chui, V. Ryzhov for valuable discussions 
and the Aspen Center for Physics where part of the results has been obtained. This work was supported by 
the National Science Foundation, Crant No. DMR 99-71518 and by NASA, Grant No. NAG8-1427. 

A. Zero sound in electron liquid 

In 1957 Landau predicted existence of collective sound excitations in Fermi systems which exist in 
coUisionless regime and called them "zero sound" (Landau 1957). Zero sound waves propagate with the 
speed Us greater than the Fermi velocity vp which excludes possibility of coUisionless Landau damping. 
Later, zero sound was observed in liquid ^He (Keen et al. 1965; Abel et al. 1966) and in electron liquid 
in metals (Burma et al. 1989; Bezuglyi et al. 1991; Bezuglyi et al. 1995). The appearance of zero-sound 
oscillations in charged Fermi liquids was studied by Silin (1958) and Gor'kov & Dzyaloshinskii (1963). In 
particular, Gor'kov & Dzyaloshinskii (1963) have shown that for weakly interacting electrons ordinary zero 
sound can be propagated in the directions of the symmetry axes in metals while spin zero sound apparently 
exists in any type of metal. As a consequence of collective motion the attenuation of zero sound can be 
atypically small compared to single particle excitations in Fermi liquids. E.g., in normal ^He liquid under 
the pressure 27atm the attenuation length of longitudinal zero sound 2500 times exceeds the free path of single 
quasiparticles and the ratio substantially increases with increasing the pressure (Corruccini et al. 1969). For 
huge pressures, such as in a neutron star crust, it is reasonably to expect very small attenuation of collective 
modes. 

Due to the presence of long-range Coulomb interactions, the Landau theory of Fermi liquid is not 
directly applicable to electron gas in metals. However, in electron gases the Coulomb interaction between 
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quasiparticlcs is screened by the gas polarizability and effectively becomes a short-range. It is therefore 
conceivable that density waves of quasiparticlcs could arise having a typical sound spectrum. If we take into 
account screening the effective interaction potential between quasiparticlcs in momentum space is given by 

= TO)' (^^' 

where £;(w,k) is the longitudinal dielectric function. In random phase approximation (RPA) (Pines & 
Nozieres 1989) 

e;(a;,k) = l-y(fc)x°(a;,k), (A2) 

where V{k) = Ane'^/k'^ is the interaction potential between bare electrons and x"(ti),k) is the density-density 
response function for noninteracting fermions which is in the limit /c -C fci?, ?jw -C Ep given by 

xV,k) = ^r^lnr^^VlV (A3) 



2Ef \2vFk \u; — vpk 

n is the electron concentration, vp is the Fermi velocity. In this section we assume absence of external 
magnetic field. Hence, in RPA 

VM^,k) = - (A4) 

[fc2 + #^(l-fln(f±i)) 

where s = ui/vpk, q^p = GTrne^/Ep is the Thomas-Fermi screening wavevector. Such form of the interaction 

potential suggests screening of the electric charge. E.g., if Vcs{k) = Aire'^ / (fc^ + q^) the interaction potential 
in coordinate space has the form Ves{r) = (e^/r) exp(— gc?'), that is the effective quasiparticle interaction is 
short-range. 

The RPA takes into account electron interaction only to the extent required to produce the screening 
field. The effects arising from interaction between screened charges (such as zero sound), are, therefore, 
beyond the RPA. One should note, however, that screening depends substantially on ui and k. For ui 3> vpk 
(s ^ 1) the screening is small 14ff(a;,k) « A-Ke^/ik"^ — q^p/2>s'^) ~ V{k) and, hence, RPA is a valid theory. 
This is the reason why RPA is successful in describing plasmons in metals for which oj ^ vpk. Indeed, the 
dispersion equation ei{uj, k) = with e;(a;, k) given by (A2) has a solution lo^ = A-Ke^n/me corresponding to 
oscillations with plasma frequency. 

However, in the region where zero sound branch is expected (for weakly interacting fermions this region 
is s — 1 <C 1) Eq. (A4) predicts large screening and RPA can not capture collective modes in this region. 
Instead, one can apply a phenomenological Landau theory based on the effective screened interactions be- 
tween quasiparticlcs. The Landau theory describes interaction by means of dimensionless phenomenological 
parameters Fi , Gi which measure the interaction strength as compared to the kinetic energy. The parameters 
Fi , Gi describe respectively the spin symmetric and spin antisymmetric parts of the quasiparticle interaction. 
For weakly interacting Fermi gas only the first two parameters Fq and Go are relevant, they are related to 
each other by Gq = —Fq which is the property of exchange interaction between electrons. Parameters Fq and 
Go describe propagation of axisymmetric (longitudinal) zero sound. Propagation of zero sound is possible 
when Fo > or Go > 0. The former case corresponds to an ordinary zero sound, while the latter one is 
known as spin zero sound in which oscillation of the quasiparticle distribution function depends on the spin 
direction. Loosely speaking Fq describes propagation through the medium of a particle-hole pair with the 
total spin 0, while in the spin zero sound the total spin of the propagating particle- hole pair is equal to 1 
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(Gottfried & Pieman 1960). The speed of ordinary zero sound Ug = svp is determined from the equation 
(Landau 1957) 



which for weakly interacting fermions Fq <^ 1 reduces to 



s- 1 « 2exp ^-2- J-^ . (A6) 

That is the speed of zero sound is exponentially close to Vp. The speed of spin zero sound is described by the 
same Eqs. (A5), (A6) with Go instead of Fq. The property Gq = — Fq suggests that one of these parameters 
is positive. This is a manifestation of a general statement that ordinary or spin longitudinal zero sound exist 
in any stable Fermi liquid (Mermin 1967). 

For a short-range weak interaction between fermions the Landau parameters can be expressed in terms of 
the interaction potential as Fq = 3nVes{k = 0)/2Ef. Formal application of Eq. (A4) at fc 0, s ^ 1 yields 
Fo ~ 2/ ln[(,s — l)/2] < 0, |Fo| ^ 1. This gives us a hint that Fq has a small negative value. The magnitude 
of Fq can be approximately estimated as the ratio of the Coulomb energy to the kinetic (Fermi) energy. For 
relativistic electron gas in the stellar crust we obtain: Fq « —e^kp/Ep = — 1/a = —0.007, where a = ch/e^ 
is the fine structure constant. Negative value of Fq suggests that Go is positive (Go = —Fq w 0.007) and, 
hence, the spin zero sound can propagate through the stellar crust. 

Such conclusion can be justified by consideration of the density-density response function, for weakly 
interacting Fermi gas it is given by (Gottfried & Pieman 1960) 

where A = 0, 1 it the total spin of the particle-hole pair which can be excited in a singlet (A = 0) or a triplet 
state (A = 1). Collective modes of the system are determined by the poles of xa('^) k). If we take Ves{uj, k) 
from Eq. (A4) then XA('^ik) has a pole at w = vpk for A = 1. This pole describes spin zero sound, while 
the plasmon branch corresponds to a pole for A = 0. In such approach we find that speed of zero sound 
is exactly equal to vp which corresponds to the limiting case of noninteracting quasiparticles. The result 
formally appears because 14ff(w, k) vanishes at s = 1 which means infinite screening of Coulomb interaction. 
However, in real systems there is no screening at distances less than the interparticle spacing « l/kp, that 
is instead of zero one should take \Vcs\ ~ e^kp/n. As a result, the position of the pole is shifted from 5 = 1 
and determined by the same Landau expression (A5) with Gq ~ e^kp/Ep. In Appendix B we discuss zero 
sound from the viewpoint of two component (particle-hole) picture of excitations in electron liquid. 

A remarkable feature of zero sound in a weakly interacting Fermi gas is that its attenuation is expo- 
nentially small. Attenuation of zero somid occurs by means of excitation of two electron- hole pairs, since 
single particle processes (Landau damping) are not allowed {us > vp). Gottfried and Pieman (1960) have 
solved the problem about zero sound attenuation in a weakly interacting Fermi gas for T = and obtained 
the following result for the imaginary part of the sound frequency 

Ima;=(s-1)^^. (A8) 



For weak interactions the multiple s — 1 is exponentially small which provides exponentially small sound 
attenuation. Formal reason for this is the logarithmic singularity in the density-density response function 
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for noninteracting fermions x°(a;,k) at uj ~ vpk (see Eq. (A3)). Finite temperature can be immediately 
included using the Landau formula (Landau 1957), finally we obtain the free path length for zero sound 

^ _us_ ^ nvpEp exp (2 + 2/ Go) 

For Go = 0.007 the exponential factor gives a multiple 10^^'', hence zero soimd propagates without attenua- 
tion across the dense stellar crust for any imaginable values of frequency and temperature at which the Fermi 
liquid theory is still valid {hu, ksT <C Ep). However, close to the stellar surface the electron gas becomes 
non relativistic. In this region the parameter Go depends on the density Go ~ e^kp/Ep = 
where as = Ti^ /rrif.e^ ~ G.53A is the Bohr radius, and could be of the order of 1 near the surface. However, 
in the surface region the effect of magnetic field becomes substantial and prevents the wave damping. We 
discuss this in another section. 



B. Equations of magnetic hydrodynamics for electron motion 

In this section we consider equations which describe wave propagation across the stellar crust. Typical 
temperature in the crust of middle-aged NSs is less then the melting temperature (Potekhin et al. 1997). 
As a result the matter forms a Coulomb crystal. One should note that zero-point ion vibrations do not 
destroy the lattice made of heavy nuclei. The amplitude of zero-point vibrations is commonly much smaller 
than typical inter-ion spacing. With increasing density the relative amplitude becomes larger, so that the 
vibrations can prevent crystallization at high p. The crystallization is completely suppressed in the so called 
quantum liquids, which exist at pe > 0.006^^Z^ (Jones & Ceperley 1996). The suppression is pronounced 
for light nuclei, however it does not occur in the lattice composed from heavy elements which is likely the 
case in the stellar crust. E.g., for iron lattice {A = 56, Z = 26) the density at which zero-point vibrations 
melt the lattice is p > 2 x lO^^g/cm^. This value is much larger than the density at the crust bottom 
(lO^^g/cm^) and, therefore, vibrations do not destroy the lattice in the crust. 

In our model we assume that crust consists of nuclei which form a lattice and a gas of free electrons. 
Free neutrons appear in the inner crust. Frequencies of radiation we are interested in are much larger than 
frequency of electron-electron and electron-phonon collisions. Hence, the electron gas can be treated as 
collisionless. As we have seen in the Appendix A, in collisionless regime zero sound is a possible collective 
excitation in electron liquid without magnetic field. It is known that in a strong magnetic field, such that 
UA ^ vp (here ua = H/ ^Airp^ is the Alfven velocity for electrons, p^. is the electron density), the collective 
excitations of electrons in compensated metals reduce to magnetoplasma waves (Buchsbaum & Gait 1961; 
Kaner & Skobov 1963). Their nature is similar to magnetosonic and Alfven waves in magnetic hydrodynam- 
ics. The dispersion relations of magnetoplasma waves can be obtained from the dielectric permittivity of 
collisionless plasma in a magnetic field. Indeed, if only electrons contribute to the permittivity then (Lifshitz 
& Pitaevskii 1979) 



a;2 



£± = 1- ,., , (Bl) 



He 



where = ^/^ire^ne/me is the electron plasma frequency and ujHe = eH/rrieC is the electron cyclotron 
frequency. For this e±_ the dispersion equation for electromagnetic waves at cj ^ lohc has two solutions 
UJ = kuA/ \/'^ + u\/c^ and to = kuA cos d / y/l + v^Jc? describing the fast magnetosonic and Alfven waves 
respectively (here is the angle between the wavevector k and H). 
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In magnetic hydrodynamics the fast magnetosonic wave reduces to usual sound in the hmit of small 
magnetic field (when ua <ti Ug). Here we demonstrate that such situation remains also valid in the coUisionless 
regime with zero sound playing the role of usual sound. For this purpose it is convenient to use a definition 
of elementary excitations in electron liquid in terms of electron-hole pairs. In such picture the excitation 
with p > pp (particle) has the energy e — vf{p — Pf), while at p < pp (hole) the energy is e = vp{pp — p), 
where is the Fermi momentum. In such definition the excitation energy is positive and gives the energy 
excess above the ground state. Elementary excitations appear in pairs so that the total number of particles 
is always equal to the total number of holes. The dynamics of excitations is described by the coUisionless 
kinetic equation for the quasiparticle distribution function n(t, r, p): 



dn dn dn f 



(9n(p') dV 



+ e 



E+- (v X H) 

c 



dn Q 
9p ' 



(B2) 



dv {2-Knf 

where v = ds/dp and p are the quasiparticle velocity and momentum. The quasiparticle charge and the 

mass are given by e = |e|sign(pi? — p) and m = meSign(p — p^?) respectively. The function /(p, p') describes 
interactions between quasiparticles and determines the velocity of zero sound. For simplicity we assume that 
/(p,p') does not contain a dependence on the quasiparticle spins. The equilibrium distribution function is 
no = l/[exp{e/kBT) + 1], so that the equilibrium quasiparticle number depends on temperature. 

Let us first consider the case of zero external magnetic field. In general, the particle-particle interaction 

differs from the interaction between a particle and a hole. As a result, it is convenient to treat our system as 
composed of two components. Then dynamics of the particle (ni) and the hole (712) distribution functions 
is described by the following equations 



dni 


dni 
dr 


dni 


'at 


dp 


dn2 


_^^9n2 
dr 


dn2 


~dt 





, dni{p') , dn2ip') 

/ll(P,P) +/l2(P,P) 

f ^ f _,-, 9ni{p') 

/22(P,P) +/2l(P,P) 



(pp' 

(27r?i)3 
(Pp' 



(B3) 
(B4) 



One can seek a solution of these equations in the form ric = n^a+Sna, where 5na — {dnoa/ de)i^a{p) exp[i(kr— 
cot)] is a small correction, a = 1,2. Here the function depends only on the direction of p since dno/de 
is nonnegligible only when p « pp. For simplicity we assume /^/^(p, p') to be constant. Then using 
dno/dp = vdno/ds and d'^p = ppdQdp we obtain equations for i^a 



-ilelEv = 0, 



(a;-kv)i/i(p)-kv j [Fiiuiip') + Fi2i^2{p')] 
{oj - kv)i/2(p) - [F22MP') + F21MP')] ^+«|e|Ev = 0, 



dn' 
An 

dn' 



(B5) 



(B6) 



where Fa0 = —Anp^ J fajsidnop/ de)dp/ {2'K?i)^ are dimensionless Landau parameters and dn' means inte- 
gration over directions of p'. Eqs. (B5), (B6) have to be solved simultaneously with the Maxwell equations 
for the electric field E. However, as we have mentioned in Appendix A, in the zero sound regime the electric 
field of quasiparticles is screened by the rest part of electrons which fill the Fermi sphere and, hence, in 
Eqs. (B5), (B6) one can put E = 0. The screening eff'ect is included into Landau parameters One 
can consider the quasiparticles as moving in a medium with a very large dielectric constant e(a;,k), then 
the Maxwell equation div(£E) = Anp results in E • k = which justifies the omitting E. In our case of the 
two component particle-hole system Fn = F22, -F12 = ^21- Under such conditions Eqs. (B5), (B6) have the 
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following solutions describing longitudinal zero sound 

CcosO , Ccos9 

s — cos s + cos a 

where s = uj/kvp = Ug/vp and 9 is the angle between the wavevector k and p. The equation for s formally 
coincides with the expression obtained by Landau for a neutral Fermi liquid 

flnfi±lVl= ' (B8) 
2 \s-l) FuTFi2 ^ ' 

in which the Fermi liquid interaction parameter F — Fn =F^i2 (Dunin & Fetisov 1972). For the solution with 
the upper sign the particles and holes move with the same average velocity and the total current is equal to 
zero. Such solution also exists if there is no screening since the particle-hole motion does not produce electric 
field. Unscreened Coulomb interaction between particles is repulsive, that is Fn > 0, while particle-hole 
interaction is attractive F12 < 0. This suggests that F = Fn — F12 > and, hence, the zero sound branch 
always exists in the system (not subject to Landau damping). For the solution with the lower sign the 
particles and holes move in opposite direction producing quasiparticle current and density oscillation. Such 
solution is allowed only when there is screening of the electric field. 

Now let us discuss the effect of magnetic field. Here it is convenient to consider each of the components 
separately. For definiteness we derive equations of magnetic hydrodynamics in terms of (quasi)particle 
motion and use Eq. (B3) with the additional magnetic term. The average quasiparticle concentration Uq, 
the quasiparticle current density jg and the average velocity Vg are given by 

"^•1 = J riiSp, ■'^ ^ / ^'^'^'^'^^P' ^"^ ^fn^d^p ' ^^^^ 
Integrating Eq. (B3) over d'^p we obtain the continuity equation 

^+VKVg) = 0. (BIO) 

To derive an equation of momentum conservation in collisionless regime one should make an assumption 
about deformation in the distribution function. The assumption is necessary only when we integrate the term 
vdni/dr in the kinetic equation. Such a term is nonnegligible as compared to the magnetic field contribution 
only in the zero sound regime. Hence, without loss of generality one can assume that perturbation in the 
distribution function is similar to those for the longitudinal zero sound. Deviation of the sound speed Ug 
from vp is determined by interactions between electrons /. However, at large densities the interaction energy 
is small compared with the electron kinetic energy and, therefore, Us ~ vp, that is s — 1 <C 1. Under such 
condition, according to Eq. (B7), the deviation of the distribution function 5ni from its equilibrium shape 
is extremely anisotropic. As a result, one may assume that drii is nonzero only in a small vicinity of the 
direction given by the vector dni/dr. Also one can omit the small terms with / describing interactions 
and take no instead of ni in the term with the electric field E. Further, we assume that w <C lvhc- Such 
assumption makes possible to omit the term with electric field in the kinetic equation. The point is that the 
magnetic field H does not allow the electron to move in a direction perpendicular to H and change its average 
momentum in this direction under the influence of E. So, if we are interested in the electron motion on a 
time scale larger than the Larmor period (that is w <C to He) only component of E along H can contribute 
to the change in electron momentum. Formally this follows from an equality E(9no/9p) = Ev{dno/de) in 
which the electron velocity v =ds/dp after averaging over fast Larmor precession has nonzero component 
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only along H. Therefore, in the kinetic equation one can substitute E by its component along the magnetic 
field. For simplicity of derivations we omit completely the term with the electric field bearing in mind 
that solutions of final equations are valid when they posses a property E • H = 0, which is the case for 
magnetosonic waves we are focusing on. Then multiplying Eq. (B3) on v and integrating over d^p we obtain 
the following analog of Euler's equation in coUisionless regime 

+ 4Vn, - — (j,xH) = 0. (Bll) 

Further wc note that if quasiparticlcs move with an average velocity then such motion induces an electric 
field E = — (l/c)(VgXH). The reason is the large electric conductivity a of the system which requires the 
total force exerting on quasiparticles to be equal to zero. The induced electric field also acts on the rest part 
of electrons filling the Fermi sphere. As a result, these electrons must move with some average velocity V to 
insure the equality E+(l/c)(V x H) w 0. Wc obtain that magnetic field couples the motion of quasiparticles 
and electrons inside the Fermi surface. Equation ( x H) ft! ( V x H) suggests that the quasiparticle current 
jg is related to the entire current density j by (jgXH) = (ng/ne)(j x H), where is the total concentration 
of electrons in the system. 

Further, using Maxwell's equations curlH =47rj/c + (l/c)9E/9t, curlE = — (l/c)9H/9t and j —a\E + 
(l/c)(VqXH)] we obtain equations of magnetic hydrodynamics for quasiparticle motion (in linearized form): 



- = curl(V, X H) + — 



— = curl(V, X H) + — ( AH--^ ) , (B12) 



^ ^ - J V". - J^(H X c„M, - ^ („ . (fg. X h)) . (^^ . h) , ,B.) 

where Pe = fneUe is the total density of electrons. The last term in Eq. (B12) describes attenuation of 
magnetosonic waves due to finite conductivity a. It is interesting to note, that for waves propagating with 
velocity close to the speed of light (which is the case for ultrarelativistic electrons in the stellar crust) the 
dissipative term goes to zero and attenuation of magnetosonic waves is substantially suppressed. 

In the limit of large conductivity {a oo) Eris. (BIO), (B12), (B13) reduce to usual equations of nondis- 
sipative magnetic hydrodynamics in which speed of sound is equal to the Fermi velocity and quasiparticle 
density undergoes oscillations instead of the total density of electrons: 



5h 
'di 



— = curl(V, X H), (B14) 



dt 



' + nqdivVg = 0, (B15) 



dt Uq ' 47rpe 47rc2pe V \ dt 

where H, Uq are equilibrium values of the magnetic field and quasiparticle density, while h, are small 
perturbations. The last term in Eq. (B16) appears due to the displacement current in Maxwell's equations. 
Similar derivations can be done for the hole component, they result in the same equations (B14)-(B16) 

with rih, V;i instead of n,, V,. Hence, one can seek a solution describing motion of one of the components 
independently, the equations of motion of the other component will be automatically satisfied, e.g., by taking 
nh = Uq and V/j = V,. 



-31 - 



Let us now derive an equation for propagation of magnetosonic waves. Taking div from both sides of 
Eq. (B16), using Eqs. (B14), (B15) and identities 

div(H X curlh) = curlh • curlH — H-curl(curlh) Ri H-Ah, 
a.(H.(^.H)).-H..„(i^.H).-H.^, 

we obtain 

The magnetic field H creates the space anisotropy in our problem. It is convenient to decompose h and Vg 
into longitudinal (parallel to H) and transverse (perpendicular to H) components: 

h = h||+hi, V, =V||+Vi. 

Then from Eqs. (B14), (B16) we find 



dt 
9V|| 



-ifdivV_L, (B18) 
(H.V||)Vi, (B19) 



3t 



Using Eqs. (B15), (B18) and (B20), we obtain 



where A|| = V^. Finally, Eqs. (B17), (B21) result in the following equation for wave propagation 

where ua = H/^A-Kpe is the Alfven velocity determined by the total electron density. 

Let us consider a plane magnetosonic wave with the wavevector k: n'^ oc ex.p{iu}t — ikr). Then Eq. 
(B22) gives the following expression for the speed of fast magnetosonic wave 



4 + + — 



J cos^ ^ ^ "'^ 



-4(1 + ^) v%u\ cos2 I /2(1 + u'Jc'), 



(B23) 

where is the angle between the wavevector k and the magnetic field H. Near the crust bottom ua ^ vp 
and u K vp (the wave reduces to zero sound) while at the star surface ua 3> vf and u « uaI \/\ + u\/c^ 
(magnetoplasma wave) . We conclude that everywhere inside the stellar crust the speed of fast magnetosonic 
wave is approximately equal to the speed of light and, therefore, trajectories of rays are approximately 
straight lines. 
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C. Propagation of electromagnetic waves across the stellar atmosphere 

In this section we discuss conditions at which electromagnetic radiation produced by vortices propagates 
across the stellar atmosphere without attenuation. The NS atmosphere is a thin plasma layer. Due to huge 
gravitational field g ^ 10^^ — 10^'''cm/s^ at the surface, the NS atmospheres are extremely compressed. The 
scale height, ~ ZhsTg/Mg {M and Z is the mass and charge of ions, Tg is the surface temperature), lies in 
the range 0.1 — 100cm (Pavlov et al. 1995). Typical densities of the NS atmospheres are 0.1 — lOg/cm^. The 
density grows with increasing g and decreases with increasing the atmosphere temperature. To be specific 
we assume ionized hydrogen atmosphere {Z = 1, m = rUe, M = nip). A strong magnetic field of the order 
of 10^^ — lO^^Gs can exist near the star magnetic poles. In such magnetic field the cyclotron frequency of 
electrons is (we take in numerical estimates H = Ax lO^^Gs) 



while the ion cyclotron frequency 



UJHe=^-^ = 1.1 X IQl^s-l, (CI) 

mc 

a;/fi = ^=6xl0iV\ (C2) 

here e is the electron charge. The strong magnetic field drastically distort the structure of atoms and increases 
the binding energies of electrons in atoms. For instance, the field lO^^Gs amplifies the binding energy of 
hydrogen atom to ~ 150eV w 1.7 x lO^K. For lower temperatures the magnetized atmospheres become 
partly ionized. In this section we assume that the atmosphere temperature is greater than lO^K so that the 
magnctoactivc plasma is highly ionized with Up = Hi = n {n^ and rij are electron and ion concentrations). 
In numerical estimates we take = 2 x 10^'K. 

An effective collision frequency of particles characterizes the plasma. For parameters of the stellar 
atmosphere whs ^ kBTs/h ~ 4 x lO^^Hz that is thermal energy is much less then the spacing between 
Landau levels of electron in magnetic field. Therefore, collisions of an electron with other particles are 
substantially suppressed because the electron motion is effectively one dimensional (only momentum along 
the magnetic field can be changed during the collision). Contrary to electrons, spacing between Landau levels 
of ions is small in comparison with the temperature {huHi <^ ksTs) and, as a result, the ion-ion collisions 
are dominant. The effective ion-ion collision frequency can be estimated as (Ginzburg 1970) 

7re4 _ , / fcsT, A 3.9n, [7^ ( 220tA 



1/3 



here Vi is the average ion velocity. Magnetic field does not afi'ect the nature of collisions provided rn ^ ro, 
where rn = Vi/iDHi is the radius of curvature of a particle in the field and td ~ Vi/ujpi is the Debye radius. 
From this we arrive at condition ujpi ^ uim- If this condition is not satisfied (as in our problem), this 
fact, however, affects only the logarithmic factor in (C3) and Eq. (C3) still can be used for approximate 
estimation of the collision frequency. For ni — 10^'*cm~^, Ts ~ 2 x lO^K we obtain that VcS lies in the 
infrared band: v^s « 4.7 x lO^^Hz. For colder NSs the atmosphere temperature is less than ~ lO^K, plasma 
is weakly ionized and ion-atom collisions become dominant. 

For n = lO^^cm""^ the plasma frequencies are 



^pe 



= 9 X 10^V\ 



TO 



4776^71 



= \l = 2.1 X 10"s-\ 
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We consider electromagnetic waves in plasma under the condition oj ^ i^eS (collisionless plasma) which 
suggests their small damping. The opposite limit uj <^ VcS corresponds to magnetic hydrodynamics. We 
choose a coordinate system so that the z axis is oriented along the magnetic field H. The dielectric permit- 
tivity tensor of a cold collisionless plasma in a magnetic field has the form (Lifshitz & Pitaevskii 1979) 



Sxx = Syy = e±, ezz=s\\, (C4) 

^xy ^yx ^V: ^xz ^yz 0, (^^) 

e^ = l-—fk_ _fk^^ (C6) 

"~ Ij^ ' ^ ' 



The dispersion equation for electromagnetic waves 



,2 



2 ^ 
k 6a0 -kakfj jSocfj 



= (C9) 



has two solutions. One of them describes fast magnetosonic waves and for frequencies not too close to 
cyclotron frequencies WRi, t^He has the form 

OJ = (CIO) 



The other solution corresponds to Alfvcn branch. General expression for the dispersion relation of this 
branch reduces to a simple form only in limiting cases. E. g., in the frequency range i/eff "C a; <C u)Hi,ijJpe 
the dispersion equation for Alfven waves is 

uj= , =, (Cll) 

'J (fccsint?)^ + 



where i? is the angle between k and H. 

Further we consider the fast magnetosonic branch only because this type of waves is generated by 
vortices. Ion-ion collisions produce damping of electromagnetic waves. To estimate parameters at which the 
atmosphere is transparent for fast magnetosonic waves one can use the following expression for the complex 
dielectric permittivity (Ginzburg 1970) 



At frequencies not too close to cyclotron frequencies the real part of the refractive index is Ref « 1, while 
the imaginary part is 

If radiation propagates through the atmosphere along the z axis its amplitude decreases by a factor exp(— u; / Imtdz /c) . 
In estimates we assume an isothermal atmosphere. Deviation of the atmosphere from an isothermal one 
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does not change the result substantially. Then concentration of particles changes with the height z as 
n(z) = no eyiY>{—Mgz/kBTs), where no is the concentration at the atmosphere bottom (z = 0). The effective 
collision frequency v^s n{z), and, therefore, Ivat oc n^ also decreases with the height according to the 
exponential dependence. Using the criteria u Jlmtdz/c w 1 one can estimate the maximum density at the 
atmosphere bottom for which the atmosphere is still transparent for vortex radiation 

Comparing Eqs. (C14) and (C3), we finally obtain 

M./V/^Ty'g < 1, 

mV4^3.97rln(220Ts/nV3)cfcs I L^/oJm, w > wj/i 

In the limit (j <^ ujhi the maximum value of no is independent of w. If g = IQi^cm/s , = 2 x lO^K 
and H = Ax lO^^Gs we obtain no ~ 4.4 x lO^^cm"^. In this limit Eq. (C15) approximately coincides with 
the photosphere density obtained from numerical calculations assuming free- free opacity (Lai et al. 1992): 
np\iot{H) « 1.8 • lQ'^^g\^^T^^'^Hi2 = 4.8 x lO^^cm"^. For w > ujHi the maximum value of no becomes 
frequency dependent and increases with increasing w. We conclude that if the particle concentration at the 
atmosphere bottom is less than 4 x lO^^cm"^ (which is usually the case) the atmosphere is transparent for 
vortex radiation (apart from narrow regions near cyclotron frequencies). 
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Fig. 1. — Schematic cross section of a dense neutron star. The star is subdivided into the atmosphere, the 
outer crust, the inner crust and superfluid outer and inner cores. 
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Fig. 2. — Generation of zero sound waves by vortices at the crust-core interface. Due to drag effect the 
hehcal vortex motion causes density oscillation of electrons in the stellar core. No sound is radiated inside 
the core because the wave length of sound is larger then the helix period. However, at the crust-core interface 
the density oscillation excites spherical sound waves that propagate across the stellar crust. 
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Fig. 3. — Oscillations of electric and magnetic fields in the fast magnctosonic wave in the limit ua S> (/,s. 
Oscillating magnetic h and electric fields E are perpendicular to each other and the wavevector k. The 
electron velocity V lies in the kH plane and perpendicular to H. 
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Fig. 4. — Mechanism of a neutron star radiation. Thermaly excited helical waves of neutron vortices in the 
superfluid core produce magnetosonic waves in the stellar crust. Magnetosonic waves propagate across the 
crust and transform into electromagnetic radiation at the star surface. Mainly the radiation comes out from 
regions with strong magnetic field (near magnetic poles). 
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Fig. 5. — Broadband spectrum of PSR B0656+14 (Koptsevich et al. 2000) and the Vela pulsar (Pavlov 
et al. 2001, Nasuti et al. 1997). Solid line is the fit by the sum of the vortex and thermal components. 
Thermal radiation of the stellar surface dominates in the ultraviolet and soft X— ray bands, while the vortex 
contribution (dash line) prevails in infrared, optical and hard X— ray bands, where its spectrum has a slope 
a » —0.45. In the far infrared band the spectrum of vortex radiation changes its behavior and follows 
Planck's formula with P(w) oc w^. 
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Fig. 6. — Broadband spectrum of the Geminga pulsar. Filled circles are data of ROSAT and ASCA (Halpcrn 
& Wang 1997), while crosses combine measurements of the Keck Observatory (Martin et al. 1998) and the 
Hubble Space Telescope (Mignani et al. 1998). Solid line is the fit by the sum of the vortex (dash line) and 
thermal contributions. 
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Fig. 7. — Temperature distribution in a neutron star crust for a model with magnetic isolating layer of width 
d = 50m. The main temperature drop occurs near the crust bottom where temperature decreases from T = 
8 X lO^K (at the corc-crust interface, z = 0) to lO^K. The rest part of the crust is approximately isothermal 
apart from a thin surface layer where temperature drops due to substantial decrease of thermoconductivity 
at low densities. Near the crust bottom the neutrino emissin results in energy losses and only a small fraction 
of the original heat flux reaches the stellar surface. 



